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SYNOPSIS 


K.V. DESIICACHAR 
Ph.D. 

Department of Electrical Engineering 
Indian Institute of Te chnology, Kanpur , INDIA 

April, 1983 

ULTRA HIGH SPEED PROTECTIVE RELAYING SCHE^ffiS 
FOR EHV/UHV TRANSMISSION LINES 
BASED ON TRAVELLING WAVE PHENOMENA 


It has long been recognized that the ultra high speed 
(UHS) clearing of faults on a transmission line improves the 
transient stability. The fault clearing time being dependent 
on the speed of the protective relay as well as that of the 
associated circuit breaker, and with the emergence of UHS 
circuit breakers, the need for developing UHS protective 
relays has become imperative. The realization of UHS 
protective relays, for the protection of EHV/UHV transmission 
lines, has been possible with the utilization of the 
travelling wave characteristics. 

The development of protective relays, based on 
travelling wave phenomena, is of very recent origin. The 
first travelling wave relay was developed by ASEA, Sweden, 
and was installed on Bonneville Power Administration’s (USA) 
500 KV system in April 1976, Since then, a few other 



travelling wave relay schemes have been proposed. Almost, 
all of these schemes necessitate the use of a fast communi- 
cation channel for internal fault detection. And, two 
of these schemes possess the fault locating feature, but 
involve in cumbersome computational procedures. Therefore, 
there is a need for a travelling wave relay scheme, v^'hich 
can detect an internal fault without the aid of a carrier 
communication channel so that the operating time of the 
carrier communication equipment is eliminated, and which 
incorporates a simple fault locating feature. 

For the development of high speed and UHS protective 
relay schemes, an accurate determination of the complex 
post-fault current and voltage waveforms is essential. For 
this purpose, accurate modelling of power systems and fault 
analysis techniques, suitable for digital simulation, are 
necessary. So far, only Johns and Aggarwal have developed 
techniques for the analysis of faulted power systems 
incorporating exact models of -transmission lines. However, 
these techniques would be applicable, with ease, only to 
simple power systems. 

It is well-known that the use of digital computers, 
for protective relaying purposes, offers several advantages. 
The digital distance algorithms, proposed for the protection 
of transmission lines, require the elimination of d.c. offset 
and also a wide spectrum of non-fundamental frequency 



components with the help of digital filtering and data 
fitting techniques, a task that entails considerable time 
delay if adequate accuracy were to be achieved. On the 
other hand, the adoption of travelling wave techniques, 
for the digital protection of transmission linos, needs the 
elimination of only the profault components, since the 
travelling wave relays use only the fault-generated compo- 
nents, Also, only a short data window is required, and 
this loads to a higher speed of fault detection. Very 
few papers have dealt with the development of digital 
travelling wave relay algorithms. Takagi et al, have 
proposed a digital travelling wave relay, which involves 
in the exchange of quantitative information between the 
ends of the protected line and consequently, imposes 
demanding requirements on the ccanmunication channel. Hence, 
there is a need for developing simpler digital travelling 
wave relay algorithms. 

Accordingly, the primary objectives of this thesis 
have been as given below. 

1. To generate a simple mathematical model for a long 

EHV/UHV transmission lino and to develop fault analysis 
techniques, which are suitable for frequency domain 
methods of computing transients in complex multinode 
power systems. 



2* To propose a new and simple travelling vjave relaying 
scheme for the protection of EHV/UHV transmission 
lines, 

3. To develop a new travelling wave relaying scheme, 
for the protection of EHVA^HV transmission lines, 
which can detect internal faults without the aid of 
a coramuni cation channel and v;hich embodies a simple 
fault locating technique, and 

4, To develop simple digital relaying algorithms, based 
on travelling wave phenomena* 

An outline of the work reported in the thesis is given 
below. 

Chapter 1 presents a brief and critical review of the 
important literature pertaining to the evolution of distance, 
travelling-wave and computer relaying schemes developed for 
the protection of transmission lines. 

Chapter 2 starts with an overview of some important 
aspects of generating the distance relay characteristics, 
suitable for tho protection of long and heavily-loaded EHV/UHV 
transmission lines, vjhich are given in the literature, Tho 
theory, principles of operation, details of hardware implemen- 
tation and test results of two new distance relaying schemes, 
fabricated with Integrated Circuits and capable of generating 
a suitable threshold characteristic, are also presented in 
this chapter. 



In Chapter 3, an accurate frequency-domain Pi model 
of a transmission line and frequency-domain nodal analysis 
techniques, suitable for the analysis of faulted multinode 
power systems, are developed. An important feature of 
the fault analysis techniques is the use of frequency- 
domain bus admittance matrix. These form an important basis 
for the digital simulation of faulted EHV/UHV networks. The 
utility of the Pi model of the line and the application of 
the fault analysis techniques are illustrated by computing 
the voltages and currents at a chosen point in tv^o sample 
power systems for symmetrical three-phase and also for one 
phase to ground faults. 

Chapter 4 presents the development of, and test 
results pertaining to, two nev; travelling wave relay 
schemes. In one relay scheme, the amplitude comparison of 
each of the three pairs of modal relay input signals is 
utilized to distinguish betv^een the reverse and forward 
faults. Tripping is initiated only if the fault is found 
to be forward at both the ends of the protected line. 

In the other scheme, one relay input signal (for each 
mode) is used for detecting reverse faults. Under reverse 
fault conditions, tripping is blocked at the remote end, 
with the aid of a fast communication channel, and also at 
the local end directly. By making use of the first, and if 
necessary the second, time derivatives of the other input 



signal, (for each mode), the instants of the first and 
second incidences of the backv/ard travelling waves, at the 
relaying point, are determined. The time interval betv^/een 
these two instants is equal to twice the travel time of the 
waves between the fault and relaying points. Thus» the 
detection of an internal fault as well as its location are 
accomplished without the aid of exchanging information between 
the two ends of the protected line. 

Both the relaying schemes cater to all types of faults 
and hence, can bo regarded as a novel type of polyphase 
relays. The underlying principles of both the schemes- are 
validated by the results obtained through digital simulation 
studies conducted on sample pov>/er systems for various condi- 
tions and two types of faults. 

In Chapter 5, the algorithms for the digital computer 
application of the amplitude-comparison and fault-locating 
relay schemes, described in Chapter 4, are developed. The 
cyclo-to-cyclo comparison method, proposed by Mann and 
Morrison, is adapted as a simple digital filter for removing 
the prefault components. The viability of these algorithms 
is tested on the digital computer, DEC System 1090 at 
I.I.T., Kanpur, by using realistic fault data obtained from 
the digital simulation of sample power systems. The 
results of those tests have demonstrated the viability of the 



proposed algorithms* In the beginning of the chapter, 
a brief overview of digital distance algorithms is 
presented in order to bring out their relative merits 
and demerits. 

The thesis concludes with Chapter 6,- wherein a 
brief review of the work carried out in this thesis and 
suggestions for further v;ork are given. 



CHAPTER 1 


INTRODUCTION 


1,1 MOTIVATION 

i 

Tl'ie ever increasing size and interconnection of power 
systems, with resultant complexity, imposes a requirement 
for protective schemes that are not only highly reliable, 
sensitive and selective, but extremely fast as well. Being 
an important link, especially in interconnected systems, a 
transmission line needs special attention with regards to 
protection requirements* The sensitivity and selectivity 
of a protective scheme depend upon the type of the relay 
unit employed. Of the several electromagnetic measuring 
units, the induction cup unit proved to be the best in 
distance relaying applications inasmuch as it was faster 
(3-5 cycles) and more sensitive compared to the other types 
of electromechanical relays and also as it was easy to 
produce any type of conventional threshold characteristics. 
However, with the advent of solid state devices such as 
semiconductor diodes and transistors, a trend towards 
employing them for relays has emerged. The need for a 
faster raeasuring unit gave impetus to the development of 
solid state relays in the initial stages. The amplitude 
and phase comparison principles, developed in its wake, gave 
rise to much greater flexibility with which special threshold 
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characteristics, like conic and quadrilateral, could be 
produced. In addition, the static relays possess other 
valuable merits such as greater sensitivity, lower burden, 
absence of contacts and mechanical motion, and immunity 
from vibrations and shocks due to external causes. The 
static relays are being used increasingly in recent years, 
particularly on EHV/UHV systems where increased sensitivity, 
reliability and speed are inportant. Recent advances, made 
in tile Integrated Circuit (IC) technology, have resulted in 
the availability of cheap and reliable IC chips. Use of ICs, 
in place of transistors and several discrete components, 
leads to compact, simpler, more reliable and economical 
relay units. Research and development work, in the field 
of relays employing ICs, is in progress. 

The selectivity, provided by a protective relaying 
scheme, depends upon the type of threshold characteristics 
obtained from the relay units employed in it. The 
selectivity, between the internal and external faults, can be 
achieved by -die use of multi- zone directional distance 
relaying schemes vjiih or without carrier-current pilot 
schemes. The selectivity, between internal faults and other 
abnormal conditions such as power swings etc., depends upon 
the shape of the threshold characteristic. So does the 
change of balance point due to fault resistance. The 
quadrilateral characteristic has piroven to be the best in 
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fulfilling these selectivity requirements to the maximum 
extent- All these factors gave rise to the motivation for 
developing new distance relaying schemes, which employ ICs 
and yield a three-step quadrilateral pickup characteristic* 

It has long been recognized that ultra high speed 
clearance of faults on interties in large power systems 
very effectively improves the transient stability limit 
(i*e* the power transfer capability for a given stability 
limit). The rotational kinetic energy introduced into a 
power system during a fault is proportional to the square 
of the fault clearance time. Therefore, high speed 
clearance of faults close to large sources of generation will 
reduce the system acceleration more than any other form of 
dynamic control which can be used only after the system is 
being accelerated. Also, it has been reported that signifi- 
cant savings in transmission investment would accrue from 
the application of high speed circuit breakers and relays 
capable of clearing faults in one-cycle primary time [l]. 

In the recent years, this aspect of improving transient 
stability has been drawing the attention of quite a few 
research engineers and organisations. The fault clearance 
time depends on the speed of the protective relay as well 
as on that of the associated circuit breaker. With the 
emergence of ultra high speed circuit breakers [l,2], the 
realization of ultra high speed protective relaying schemes 
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has become imperative* 

An accurate computation of the relaying quantities 
during a very short interyal of time immediately following 
fault occurrence is essential, in order to design and assess 
the perfoiroance of not only analog ultra high speed protective 
schemes but also digital computer relaying schemes* Modelling 
of the power system for this purpose needs consideration of 
the distributed nature of the transmission line parameters, 
the nature of the earth return, the frequency dependency of 
the line and earth parameters, and the effect of untransposed 
lines, on account of the fact that high frequency electro- 
magnetic transients, in the form of travelling waves, exist in 
the post fault waveforms during this short interval of time* 
The frequency domain method of computing these electro- 
magnetic transients, although cumbersome, is the most 
accurate one* So far, as it appears from the recent 
literature, only one paper [3] has dealt with the frequency 
domain analysis of power systems, with all these factors 
taken into consideration. However, the mathematical formula- 
tion, developed therein, is suitable only for the analysis of 
simple power systems. From this, it is evident that there is 
a need for a more generalized mathematical formulation, which 
would be, suitable for the analysis of complex multinode power 
systans* ^ 
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The quest for an ultra high speed protective scheme, 
which is unaffected by the electromagnetic transients, has 
culminated in the generation of a novel and strikingly 
simple relaying concept* a concept that makes use of the 
characteristics of these very electromagnetic transients, 
manifested as travelling waves, to distinguish faulted condi- 
tions from noMial as well as all other abnormal conditions, 
and internal faults from external faults. Only a few ultra 
high speed relaying schemes, embodying the above concept, 
have been proposed so far. However, there is a need to 
search for better alternative schemes based on this concept 
and to examine their viability in order to provide the 
practising engineers with a wider, if not better, choice. 

Digital computers have been applied to, or considered 
for, multifarious off-line and on-line tasks pertaining to 
power system analysis, control and operation. Their use for 
the protection of power system equipment is of recent origin, 
the first serious proposals appearing in late 1960's. Also, 
very recently, there has been a trend towards employing 
microprocessors and multiprocessors for power system prote- 
ction purposes. Uptill now, quite a good number of computer- 
oriented real-time distance protection algorithms have been 
developed. Extraction of the fundamental components from 
the complex post fault waveforms with iiie aid of analog and/or 
digital filters is necessary with the above algorithms. But, 
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use of these filters entails time delay. Contrarily, the 
algorithms based on travelling wave relaying concept^^ 
need no elaborate filtering and consequently, they would 
be both faster and accurate. It appears that only two 
computer relaying schemes, employing an algorithm of the 
latter type, have been proposed so far. Therefore, there 
is certainly a need for better and simpler algorithms of 
this type. 

1.2 OBJECTIVES AMD SCOPE 

The objectives and scope of the work reported in -Uiis 
thesis have been j 

1. To present a critical review of the important solid- 
state relaying schemes reported so far, to explore the 
ability of the solid-state circuitry, in terms of simplicity, 
reliability and flexibility, to generate any desired complex 
pickup characteristics and to develop a faster and more 
reliable digital relay which uses digital ICs, and which 
generates a multizone quadrilateral pickup characteristic, 
acknowledged to be the best for the protection of long and 
heavily loaded EHVA^HV transmission lines. 

2. To develop a simple mathematical model for a long 
EHV/UHV transmission line and techniques of fault analysis 
based on nodal formulation which are suitable for frequency 
domain methods of computing transients in complex multinode 
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power systems. 

3. To present an overview of tlie digital computer relay- 
ing algorithms developed uptill now as well as the philosophy 
behind the existing and proposed new algorithms and 
schemes. 

4. To develop new relaying principles based on travelling 
wave phenomena, and computer relaying algorithms based on 
these relaying principles; to develop a travelling wave 
relay which incorporates the fault locating capability and 
which dispenses with the need for carrier or microwave 
communication channels in the event of internal faults so 
that the fault detection time is reduced. 

1.3 LITERATURE SURVEY 

In the early days of electric power transmission, over- 
current relays were mostly used for the protection of 
transmission lines. However, as higher and higher trans- 
mission voltages, necessitated by the increasing size of the 
power systems, were employed, these o.c. relays were found 
unsuitable because of the following reasons. The balance 
point of an overcurrent relay shifts with the type of fault, 
and also with variation in generating capacity. There may 
be a wrong tripping in the event of a fault on a parallel 
line under certain conditions [4]. Also, they can be 
used only on systems where the minimum fault current 
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exceeds the maximum load current. However, the directional 
overcurrent relays are still used for ground fault relaying 
at subtransmission and transmission levels [5]. On account 
of the above-mentioned problems associated with the overcurrent 
relays, the distance relays have been developed. In view of 
the necessity for ultra high speed fault clearance for the 
attainment of improved transient stability, protective 
relays, sNhose operation is based upon the travelling wave 
phenomena, are being developed, A brief and critical review 
of the important literature pertaining to the evolution and 
development of distance and travelling-wave relays and 
computer relaying schemes for the protection of transmission 
lines, is presented in the follov/ing sections. 

1.3.1 Distance Relays 

1.3. 1.1 Electromechanical Relays 

In the evolution of relays for the protection of trans- 
mission lines, the electromechanical relays were developed 
first, and in that of distance relays, the plain impedance 
relay was the first one that was conceived and developed. 

In 1923, Crichton [6] reported about an impedance relay which 
employed an induction disc actuating structure and operated 
in a time proportional to the impedance between the relay 
and the fault point. In 1928, McLaughlin and Erickson [7] 
reported about a directional impedance- time relay built with 
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an induction disc actuating structure. They presented 
the constructional details, and described the principle 
of operation, methods of obtaining proper voltage for the 
restraining elements and methods of compensating for the 
voltage drop in power transformers and for the secondary 
neutral shift. In 1930, Crichton [8] reported about a high 
speed impedance relay, the high speed of operation having 
been obtained with a balanced beam structure. He described 
a three-step directional impedance relaying scheme and 
discussed about the effect of fault resistance on the 
operating speed of the relay. The development of a high 
speed impedance relay with a balanced beam structure was 
also reported by Goldsborough and Lewis [9]. In 1944, 
Goldsborough [lO] reported about a modified impedance relay 
built with a balanced beam structure and described how a 
circular pickup characteristic of any desired radius and 
with any desired location of the centre could be obtained with 
a voltage proportional to the system current added to or sub- 
tracted from the relay voltage and with a current proportional 
to the system voltage added to or subtracted from the relay 
current. 

In 1928, Warrington [4] designed an induction disc type 
reactance relay for an American Company. In 1931, George [ll] 
presented the performance records of normal speed as well as 
high speed reactance relays installed at two hydroplants of 
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Tennesse Electric Power Company of America, and reported 
that, although they had performed satisfactorily on long 
lines with wide variations in short circuit conditions, 
certain difficulties had arisen on interconnectors operating 
near stability limit. In 1933, VJarrington [l2] reported 
about a high speed reactance relay which was built with a 
four-pole induction cup actuating structure carrying current 
coils on one pair of opposite poles, and current and 
voltage coils on each of other poles. Dewey and Me Glynn [l3] 
as well as Warren C. New [l4] have reported about three-step 
directional reactance relaying schemes with zones 1 and 2 
protected by the combination of reactance and mho units, and 
zone 3, by mho units. Both the reactance and mho relays 
had employed 4-pole induction cup structures. 

Though the mho relay was first used in- 1933 as the 
directional unit for an early type reactance relay [l2], its 
independent use for the protection of heavily loaded and long 
transmission lines was first recommended in 1943 by Warrington 
[l5] with its merits lucidly brought out. In 1944, Cordrey 
and Warrington [l6] described its actual use in a carrier 
current schene. Later, Hutchinson [l7] described its use in 
a three-step distance relaying scheme in which protection for 
zones 1 and 2 was provided by normal mho units, and for 
'zone 3, by an offset mho unit. He presented also the test 
results in the fom of transient overreach of zone 1 unit, 
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and the operating time of each unit as a function of the 
ratio of the impedance to fault to the setting impedance. 

In 1962, Skuderna [l8] put forth the mathematical 
development of how offset conic and limacon characteristics 
could be produced with a four-pole induction cup structure. 

He presented also detailed circuitry needed to obtain the 
various winding inputs and described the setting procedure 
with a numerical example for getting an offset elliptical 
characteristic* 

1.3. 1.2 Electronic Relays 

In 1934, V^ideroe Cl9] presented electronic circuits, 
incorporating thyratron tubes, which were equivalent of 
many electromechanical relays in use at that time. In 1948, 
Mac Pherson and Viarrington [20] described an electronic mho 
relay v>rherein instantaneous values of voltage and current 
inputs to •die relay were compared at the instant of voltage 
input maximum. Barnes and Mac Pherson [21] recounted the 
field experience of this relay and reported that it was not 
adversely affected by unknown system or other disturbances. 

In 1949, Loving [22] published electronic circuits for many 
common protective functions, and presented experimental 
results. In 1954, Bergse-th published a paper [23] describing 
an electronic directional distance relay which was insensitive 
to moderate waveform dis'tortion and which was opera -ted by 
comparing the phase angle between two derived voltages. Hodges 
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and Mac Pherson [24], Price et al« [25] and Seely and Koss 
[26] described the operating principles and performance of 
an all-electronic one-cycle carrier relaying system. However, 
with the advent of solid-state devices, the development of 
these relays has ceased. 

1, 3,1,3 Solid-State Relays 

The first serious proposal, for the employment of transis- 
torised circuits for power system protective relaying, came 
from Adamson and Wedepohl [27,28] in 1956, In the first 
paper [27], they presented the mathematical development of 
determining the inputs necessary to obtain the directional, 
ohm, offset impedance and mho characteristics with a two- 
input phase comparator, detailed circuitry of the prototype 
relays built by them, and dynamic test results for a mho 
relay. In the second [28], they presented the dual phase 
comparator technique with which the transient overreach, 
present in the coincidence-block average comparators, could 
be reduced, and substantiated its effectiveness with the 
results of dynamic tests conducted on a prototype mho relay 
incorporating this technique. In 1959, Adamson and Talkhan 
[29] described a phase-comparison carrier relaying scheme, 
that employed transistors, and presented the test results. In 
1960, Dewey and Hodges [33] described a transistorised phase - 
comparison carrier relaying scheme in which transistor cir- 
cuits were used for level detector, comparer and squaring and 
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tripping amplifiers. Dewey et al. [31] and Caleca et al.[32] 
described, in two companion papers, the design principles, 
typical circuits, application and test results of a directio- 
nal-comparison pilot relaying scheme which employed solid- 
state mho distance units with offset impedance supervision 
for phase faults and which employed the block-spike phase 
comparison technique for producing the distance relay, 
characteristics. Hahn [33] described a three-step solid- 
state directional impedance relaying scheme employing 
transistors and printed circuit boards, and the coincidence- 
block average phase comparison technique. 

In 1965, Vifedepohl [34] reported about a polarised mho 
relay which surmounted the problem of nonoperation for close- 
in faults. In the same year, Humpage and Sabberwal [35] 
presented a mathematical basis for two-input phase compara- 
tors with angular limits of comparison other than +90° and 
gave block schematic diagrams of symmetrical and unsymmetrical 
comparators, the inputs necessary for obtaining composite 
characteristics and a plural comparator arrangement that 
generates a quadrilateral characteristic. In 1966, 
Parthasarathy [36,37] reported about a static three- zone 
quadrilateral distance relay incorporating trapsistor cir- 
cuits for measuring, discriminating and logic functions, and 
employing the multi-signal block-spike phase comparison 
technique. The papers include the results of dynamic tests 
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conducted on prototype (laboratory model) relays* In the 
same year, Kans Hoel et al, [38] described how reduced~size 
third zone pickup characteristics could be produced by 
multi-signal phase comparison* 

In 1967, Anil Kumar and Parthasarathy [39] presented a 
mathematical basis for multi-input sine phase comparators and 
synthesised the inputs required to obtain conventional as 
well as quadrilateral characteristics. In the same year. 

Anil Kumar et al. [40] reported about a three-step transis- 
torised quadrilateral distance relay using the same inputs 
as in [36] and discussed the modifications necessary for using 
it with carrier-blocking and carrier-tripping schemes. In 
1968, the optimum dynamic design aspects of static relay 
comparators wore presented in a paper by Jackson et al. [41]. 

In the same year. Me Laren [42] developed a sampling technique 
which permits a comparison of instantaneous values of relay 
inputs at different instants of time and dispenses with the 
need for phase- shifting and mixing circuits. The theory and 
design of a solid-state high speed phase comparison relay 
based on the zero-crossings of the line currents were publi- 
shed in two companion papers by Robertson and Norman Meikle 
[43,44]. Vitanov [45] reported about the production of quadri- 
lateral characteristics witii transistorised circuitry. 

In 1970, Ramamoorty and Wani [46] reported about the 
fabrication and test results of a solid-state quadrilateral 
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distance relay. In 1970, Anil Kumar [47] published a 
paper describing a new sampling circuit that employed a 
combination of ferrite cores and semiconductor components 
and also a new technique that needed sampling of only one 
input, the other being converted to a d.c. level. The 
advantage of this new technique is that it does not require 
synchronisation of input signals. Khincha et al. [48] put 
forth a new amplitude comparison technique in which full- 
wave rectified, but unsmoothed, voltage and current combi- 
nations are compared at every instant during each half 
cycle, and derived the inputs necessary for generating conic 
and quadrilateral characteristics. They extended the above 
principle [49] to produce parallelogram, directional 
parallelogram and directional quadrilateral characteristics. 
Anil Kumar [50] described the production of quadrilateral 
characteristics with phase sequence detectors. He presented 
two types of phase sequence detectors} one uses ferrite core 
logic and the other, semiconductor flip-flops. 

In 1972, Johns [51] presented a generalised phase 
comparison technique in which a comparison period of fixed 
duration and related to either positive-going or negatives- 
going zero crossing is assigned to each relaying signal and 
the coincidence of these comparison periods of all the 
relaying signals provides the criterion for the relay opera- 
tion. He described how settings for two-input comparators 
could be determined. He reported [52] extension of this 
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technique to multi-input comparators and described the 
procedure to. determine the settings, with an example 
worked out for a 4-input comparator- He also described 
how a quadrilateral characteristic could be generated with a 
3-input comparator employing this technique. In i973, he 
[53] presented variable characteristic generalised techniques 
for distance protection- These techniques involve in varying 
the resistive reach of a replica impedance by adding an extra 
current^derived component to the appropriate relay input- 
This component is a linear function of all the line currents 
and is arranged to be zero under all balanced conditions and 
finite during unbalanced faults* He described how to get a 
variable quadrilateral characteristic vAxose area will be 
minimal under power swings and heavy load conditions. He 
reported also [54,55] the extension of this principle to the 
distance protection of double circuit lines and lines with 
selective-pole autoreclosure facility. 

In 1974, Jackson [56] presented a comparator with auto- 
matically variable angular limits of phase comparison. The 
automatic variation is accomplished by modifying the integra- 
tor response during the anticoincidence period with pulses of 
widths corresponding to anticoincidence periods. In this way, 
shaped pickup characteristics, narrow along the power swing 
locus, can be obtained. In the same year, Paithankar and 
Ingole [57] reported the production of quadrilateral 
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characteristics with multi-input amplitude comparators 
incorporating new techniques. In 1979, Ramamoorty and Lall 
presented a paper [58] describing a versatile phase comparator 
that was built with digital circuits and which employed 
coincidence-block average principle. Another novel feature 
of ttiis comparator is a digital up/down counter which does the 
same job as the capacitor in the integrating circuit of a 
conventional comparator. Also, asymmetrical angular limits 
of comparison can easily be obtained by having clocks of 
different frequencies during coincidence and anticoincidence 
periods. 

In 1980, Parthasarathy et al, [59] presented a new 
distance relay with an adaptive pickup characteristic which 
has a narrow tripping area during power swing conditions, 
and which automatically expands to large area during unbalanced 
faults. The adaptive feature is realised by controlling the 
angular limits of phase comparison with negative or zero 
sequence voltages generated during unbalanced faults. The 
results of the dynamic tests conducted on a poliarised mho 
relay employing this technique were also presented. In the 
same year, Ramamoorty and Lall [60] reported about a multi- 
input phase comparator using digital ICs,- A 3-input comparator, 
yielding a quadrilateral characteristic and possessing the 
capability of generating any slope for the sides of the 
characteristic, was described in the paper, A solid-state 
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distance relay, employing an operational amplifier chip as 
an amplitude comparator and producing elliptical characteri- 
stic, was developed by Ramamoorty et al, [61], Paithankar 
and Thoke [62] reported about an earth-fault quadrilateral 
distance relay, capable of producing an adaptable swivelling 
characteristic and suitable for the protection of double end 
infeed lines* In the same year, Weller et al* [63] put 
forth a new technique of phase comparison* The relay inputs 
are squarewaved, treated as binary and AND compounded. The 
states of the AND gate output will have a particular sequence 
over the operating angular range, which is entirely different 
from that over the blocking range* With the aid of a counter, 
the tripping sequence is determined. In 1982, Kellog et al, 
[64] presented a relaying scheme using digital circuits and 
CMOS ^ogic. The phase comparator described by them can be 
used as either sine or cosine comparator. It can also be 
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used, with minor modifications, as a polyphase distance relay. 

Br*4ten and Hoel [65], Ryder et al, [66] and Giraoyan [67] 
described protective relays employing rectifier bridge phase 
comparators in conjunction with a sensitive output relay, 
usually of the polarised moving iron or moving coil type* 

1.3.2 Travelling Wave Relays 

For the purpose of designing, developing and assessing the 
performance of ultra high speed relays, especially those based 
on travelling wave phenomena, fault analysis of multinode power 
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systems, with such factors as the distributed and untransposed 
nature of transmission lines and the frequency dependency of 
the line and earth parameters taken into consideration and 
suitable for digital simulation, is essential. Humpage 
et al. [68] reported a digital simulation of detailed models 
of synchronous generators, transmission lines, and power and 
protective transformers* Kothari et al. described modelling 
techniques for simulating the primary system and the trans- 
ducers in one paper [69], and digital simulation of fault 
currents and voltages in the other [70]. Johns and Aggarwal 
[3] described the digital simulation of faulted lines with all 
the factors mentioned above, taken into consideration. They 
have also described fault analysis techniques# 

Vi tins [71] described a novel transmission line protective 
scheme based on the travelling wave phenomena. Fault location 
is identified as fault detection and the location of the fault 
is determined by a time delay equal to twice the travel time of 
travelling waves between the fault and relaying points. A 
correlation technique is developed to determine this time delay* 
Yee and Esztergalyos [72], Chamia and Liberman [73], 
Esztergalyos et al* [74], Matele [75] and Garter [76] reported 
about the development, design, application and operational 
experience of an ultra high speed relaying scheme based on 
travelling wave theory* This relaying scheme utilises the 
deviations in the relaying-point voltages and currents 
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occurring as a result of a fault and compares the polarities 
of the voltage and current deviations to determine as to 
whether the fault is ahead of or behind the relaying point* 

By exchanging information between -tiie ends through a communi- 
cation channel, tripping is initiated in the event of an 
internal fault, Johns [77], and Johns and Aggarwal [78] 
described an ultra high speed relaying scheme in which the 
direction to fault is determined by the sequence in which 
relay input signals, related to the modal components of the 
relaying-point voltage and current deviations, exceed a 
preset level. This scheme is designed to operate in 
conjunction with a carrier communication channel and in a 
carrier blocking mode, 

Dommel and Michels [79] described a new method of fault 
detection which utilizes travelling wave voltage and current 
transients, A discriminant is derived and used to distin- 
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guish internal faults from external faults, Takagi et al, 
[80,81] and Akimoto et al, [82,83] presented the theory, 
sensitivity analysis and test results of a differential relay 
scheme based on travelling wave theory. They defined a 
discriminant which is a function of surge impedance of the 
line and instantaneous terminal voltages and currents. The 
discriminant is zero for normal conditions and external 
faults, and finite for internal faults, Ramamoorty and Verma 
[84] described a relaying principle based on travelling wave 
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phenomena* A discriminant^ which is a function of backward- 
wave voltage and current deviations and a setting resistance, 
is defined* If this discriminant exceeds a preset value, 
which is the case for internal faults, tripping is initiated* 
Vi tins [85] defined two wave signals, formed from the current 
and voltaga deviations, and described the determination of the 
direction to a fault by detecting which of these two signals 
first reaches a threshold constant. He presented a geometric 
approach for this purpose. 

1*3*3 Digital Computer Relaying 

The use of digital computers and microprocessors for 
protective relaying purposes has been engaging the attention 
of research and practising engineers since late 1960s. The 
first serious proposals for using digital computers came from 
Rockefeller [86]# The algorithms proposed till now can be 
classified into two categories) distance relay algorithms and 
travelling wave relay algorithms. The distance relay algori- 
thms involve in the determination of the fundamental-freQuency 
impedance to fault from the fundamental components of voltages 
and currents which are extracted from the complex post-fault 
waveforms by analogue and/or digital filters* Slemon et al. 

[87] described the determination of the fundamental components 
in phasor form from an ensemble of samples collected over one 
full cycle by using Fourier Analysis techniques* Ramamoorty 

[88] described the determination of the impedance from the peak 
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values of voltage and current and phase angle between -Oiem 
and from the fundamental phasor components of voltage and 
current, both these methods requiring samples collected 
over one full cycle* Mann and Morrison [89] described the 
predictive calculation of the peak values of and the phase 
angle between the voltage and current from a much fewer 
number of samples and their time derivatives and also dis- 
cussed, in another paper [90], the relaying of a three-phase 
line using this technique* The predictive technique was 
presented also by Gilcrest et al* [91] who proposed the use 
of first and second time derivatives of the samples, and by 
Gilbert and Shovlin [92] who proposed the use of samples and 
the sampling interval. Carr and Jackson [93] described the 
use of tvro digital orthogonal notch filters, with sine chara- 
cteristics to determine the magnitude and phase angle of the 
fundamental components from samples taken at four equally 
spaced time intervals over the fundamental period, and the 
use of an analog low pass RC filter with a cutoff frequency 
of 85 Hz to band-d.imit the voltage and current signals prior 
to sampling. 

The use of Fourier Analysis techniques for extracting the 
fundamental components in phasor form has be.en reported by 
several authors, besides Slemon et al* [87] and Ramamoorty [88]. 
Hope and Umamaheswaran [94] described the use of odd and oven 
square waves, besides sine and cosine functions, for the 
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extraction of the fundamental components. Hope et al. [95] 
and Tirupathaiah et al, [96] have also employed Fourier 
Analysis techniques with one cycle data window, whilst 
Phadke et al. [97,98] and Wiszniewski [99] have employed 
these techniques with half a cycle data window. Horton [ICX)] 
described 1iie use of Walsh functions to extract, the fundamen- 
tal components. 

Quite a few authors have proposed the modelling of the 
line by a series R-L circuit, and solving the resultant 
differential equations to obtain the values of R and L between 
the relaying and fault points. Me Ihnes and Morrison [lOl] 
proposed the integration of the differential equations over 
two successive sampling intervals to generate adequate 
number of equations to solve for R and L, and used the 
trapezoidal rule to evaluate the integrals. Poncelet [l02], 
and Bernard and Bastide [103] treated the deficiency in this 
modelling as an error and solved for R and t subject to the 
minimisation of the sum of squares of these errors over a 
certain number of successive sampling intervals. Smolinski 
[104] modelled the line by a single Pi section and solved 
the resultant differential equations by replacing them by 
finite differences and using four sample sets to obtain the 
values of R and L. Ranjbar and Cory [105] presented a novel 
digital haimonic filter in \Ahich the integration is carried 
out over a certain number of overlapping subintervals with 
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preset limits while the above differential equations are 
being solved by the numerical integration method. In this 
way, any unwanted harmonics and their multiples can be 
eliminated. Sakaguchi and Uemura [106] defined the inverse 
Laplace transform of l/(sL+R) as a weighting function and 
discussed a numerical solution technique to determine it. 

They showed that the weighting function is greater than 
zero only for a fault ahead of tbe relaying point. 

Sachdev and Baribeau [107] assumed the post-fault 
wavefoim to be comprising of a decaying d.c. offset, the 
fundamental and a certain number of harmonics and determined 
the unknown parameters of the fundamental by applying the 
least square error method. After expanding the exponential 
term associated with the d.c. offset in a series, they 
considered three terms of this series and only the third 
harmonic for elimination. Following the same method, Luckett 
et al. Cl08] considered the d.c. offset and two harmonics for 
elimination, v\tiilst Brooks [l09] represented the complex post- 
fault waveform by a constant plus the fundamental frequency 
term. Johns and Martin [llO] presented a finite transform 
method by which the fundamental frequency impedance can be 
determined by carrying out the filtering process in the 
frequency domain. Girgis and Brown [ill, 112] developed two- 
state and three-state Kalman filters to extract the fundamental 
voltage and current phasors respectively, and reported that 
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the error involved is less than 1 percent after half a 
cycle, Sanderson and V^right [113] have dealt with series 
compensated lines. They modelled the line by a series 
RLC circuit and solved the resultant differential 
equations by Integrating them over three successive sampling 
intervals. 

Vi tin's correlation method based on travelling wave 
theory [71] is suitable for computer application also. 

Takagi and Yamakosi [ll4] described the microcomputer 
application of the differential current relay [80--83] whose 
operation is based on the travelling wave theory. 

1.4 SUMMARY OF THE WORK REPORTED IN THIS THESIS 

The summary, of the work carried out and reported in 
this thesis, is presented below chapterwise. 

Chapter Two starts with an overview of some important 
aspects of generating the distance relay characteristics for 
protection of long and heavily-loaded EHV/UHV transmission 
lines. Later, two new distance relays, capable of generating 
a three-step quadrilateral pickup characteristic, are 
developed. Each proposed relay consists of a forward fault 
detector, which generates an open quadrilateral pickup cha- 
racteristic, and a zone detector, wrtiich closes the open 
quadrilateral characteristic at the appropriate zonal reach, 
depending upon the location of the fault. 
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Both the relays are fabricated with the ICs. The 
first relay scheme demonstrates that the use of ICs 
greatly simplifies the circuitry, although the input 
signals used are the same as given in reference [36] • The 
second relay has been fabricated only with the digital ICs, 
a fea-feire that enhances the accuracy of measurement* Also, 
its forward fault detector uses only three input signals, 
which results in enhanced reliability* Although, the manner, 
in which the open quadrilateral characteristic is produced in 
this scheme, is somewhat similar to that described in 
reference [60], the modifications incorporated herein render 
the measurement of coincidence period to be completed in 
exactly half a cycle as against more than half a cycle in 
their scheme. The relays have been fabricated and tested in 
the laboratory* 

In Chapter Three, the development of an equivalent Pi 
circuit for an untransposed transmission line, without 
disregarding the effect of a host of frequencies present in 
the complex post-fault waveforms, is described* This model 
makes the frequency- domain nodal analysis, of multinode power 
systems easy. Also, a, detailed fault analysis of multinode 
power system networks, based on nodal formulation and suitable 
for frequency domain methods, is presented in the same chapter* 
This analysis forms an important basis for the digital simul ac- 
tion of faulted e.h.v. networks* The utility of the Pi model 
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of the line developed, and the application of the fault 
analysis presented, are illustrated by determining the 
voltages and currents at a chosen point in two sample power 
systems for symmetrical ttiree phase and single line to 
ground faults. For this purpose, computer programmes are 
developed for determining the frequency-dependent parameters 
of the transmission line, and for. determining the time domain 
valued; from the frequency domain values through modified 
Fourier transform techniques. In this connection, it is 
shown that, the application of fast Fourier transform 
techniques in the computation of power system transients 
provides only marginal relief in terms of computational 
speed. . 

In Chapter Four, two new ultra high speed relaying 
schemes, utilizing the characteristics of the travelling 
waves generated immediately after the occurrence of a fault 
and for the protection of three-phase transmission lines, 
are developed. In one scheme, amplitude comparison of each 
of the three pairs of relay inputs, formed from the modal 
components of the fault-generated components of the relaying- 
point phase voltages and line currents, yields positive 
information as to whether the fault is ahead of or behind 
the relaying point. Exchanging information between the two 
ends of the line through a fast communication channel, 
tripping is initiated if the fault is found to be ahead of the 
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relaying points at both the ends of the line. In the second 
scheme also, two relay inputs for each mode of propagation 
are used. But, in this scheme, the magnitude of one input 
is used only to determine whether the fault is behind the 
relaying point. If the fault is behind, operation of the 
relay scheme is blocked, directly at the local end, and 
through a communication channel, at the remote end. The 
first and, if necessary the second time derivatives of the 
other input are used not only to detect an internal fault, 
but also to determine its distance from the relaying point. 

The reflection properties of the travelling waves are made 
use of, and no communication between the ends is necessary, 
in the detection and location of an internal fault. Both 
the schemes cater to all types of faults and therefore, can be 
regarded as polyphase relays. The underlying principles, of 
both the relaying schemes, are validated by the results 
obtained through the digital simulation studies conducted on 
sample power systems for various conditions of faults. 

Chapter Five commences with a brief overview of digital 
distance algorithms in order to bring out their merits and 
demerits. Next, new algorithms for the digital computer 
application of the amplitude-comparison and fault-locating 
relay schemes, described in Chapter 4, are developed. The 
cycle-to-cycle comparison method, proposed by Mann and Morrison 
C 90 ], is adapted as a simple digital filter for eliminating the 
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prefault components. The viability of these algorithms is 
tested on the digital computer, DEC System 1090 at I.I.T.,, 
Kanpur. For this purpose, realistic fault data, obtained 
from the digital simulation of sample power systems, are 
used. The results of these tests have demonstrated the 
viability of the proposed algorithms. 

Finally, the thesis concludes with Chapter 6, which 
highlights briefly the work reported in this thesis and 
also discusses the future scope of work in this field. 



CHAPTER 2 


PROTECTION SCHEMES FOR EHV/UHV TRANSMISSION LINES 
USING SOLID STATE COMPONENTS 

2.1 INTRODUCTION 

The rapid growth in the size and complexity of 
electrical power systems has given rise to increasing demands 
for reliable, faster and more discriminative protective 
schemes. The electro-mechanical relays, which have rendered 
satisfactory service over several years, have not been able 
to cope up with these demands. Consequently, there has been 
a constant search for better alternatives. First, the 
electronic relays, incorporating vacuum and gas-filled 
tube circuits, have made their appearance as a possible 
alternative. However, notwithstanding the absence of 
contacts with ensuing ease of maintenance, lower burdens 
on c.t.s and p.t.s and fast in operation even when close 
to pickup level, their development was short-lived due to 
the following reasons. They were not robust, they required 
complex wiring and their filaments imposed a burden on the 
station batteries. However, with the emergence of solid- 
state components, the solid-state relays, which provide a 
much better alternative to l^e electromechanical relays, 
could be developed. Amongst the several advantages of 
solid-state relays, ability to provide optimum operating 
speeds, flexible control of the shape of the threshold 
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characteristics and reduced panel space requirement are 
worth mentioning* Recently, the development of micro- 
electronics has made available cheap and highly reliable 
integrated circuit (IC) chips* As is well known, a single 
chip integrates several complex circuits into a tiny compo- 
nent. Thus, the use of ICs eliminates the need for wiring 
complicated circuits, improves the reliability and at the 
same time, it greatly reduces the panel space requirements. 

The development of distance relciy characteristics, 
suitable for the protection of long heavily-loaded trans- 
mission lines, has been the concern of protection engineers 
since a long time. Of the various distance threshold chara- 
cteristics developed for this purpose, the quadrilateral 
characteristic is considered to be the best since it 
encloses the fault area very compactly and consequently, 
possesses the valuable properties of least tendency to 
operate under power-swing and heavy-load conditions, and 
also greater immunity to the effects of fault resistance* 

In the light of the above, an overview of some important 
techniques and methods, developed to generate distance relay 
characteristics intended for the protection of long heavily- 
loaded e.h.v* and u.h.v. transmission lines is presented 
first, in order to throw light on their limitations, merits 
and demerits* Next, the principle of operation, hardware 
implementation and test results of two new relaying schemes 
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using ICs and capable of generating a three-step 
quadrilateral characteristic are presented* 

2.2 OVERVIEW OF DISTANCE RELAY CHARACTERISTICS INTENDED 

FOR THE PROTECTION OF EHV/UHV TRANSMISSION LINES 

One of the chief requirements of the distance relay 
characteristics, meant for the protection of long heavily 
loaded e.h.v* and u*h,v, transmission lines, has been 
that there shall be no encroachment into the tripping 
area by heavy load and also during power swing conditions. 
Basically, there have been five approaches which have aimed 
at fulfilling this requirement. They are the use of 
elliptical, composite, shaped, adaptive and quadrilateral 
characteristics. An overview of a few important methods 
and techniques developed to generate these characteristics 
is presented in this section. 

Braten and Hoel [65] employed a three-input circulating- 
current rectifier-bridge amplitude comparator to produce 
normal and offset elliptical characteristics. Ramamoorty 
et al. [6l] used phase-splitting networks to rectify the 
input signals and an operational amplifier (op amp) summing 
circuit for amplitude comparison of these signals for 
generating normal and offset elliptical characteristics with 
three input signals. The round nature of these characteri- 
stics makesthe response of the relay, to resistive faults 
near a zonal reach, poor* 
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The generation of composite pickup characteristics, 
which are composed of circular arcs and straight line 
segments and which have a tripping area closely surrounding 
the fault area, has been reported by Humpage and Sabberwal 
[35j, Hans Hoel et al, [38] and Johns [51,52], Humpage 
and Sabberwal, and Hans Hoel et al, described the generation 
of a characteristic, composed of a mho circle and two 
blinders, with three signal phase comparison,. Johns 
described the generation of a characteristic, which is a 
combination of restricted directional, restricted ohm and 
restricted mho characteristics, with three signal phase 
comparison embodying a new generalised phase comparison 
technique. With this technique, symmetrical angular limits 
other than + 90° as well as asymmetrical angular limits can 
be obtained with greater facility, A pulse of fixed dura- 
tion and displaced from the positive going (or negative going) 
zero crossing by a fixed time delay is assigned to each 
relaying signal. The coincidence of such pulses of all the 
relaying signals is measured to provide a trip signal. Then, 
the angular limits would be dependent both on the pulse 
widths and time delays. For example, in 'Uie case of a two- 
input comparator, the limits are given by 

(e^- 02 + 6 ^) where 6 ^^ and 6 ^ are, respectively, the pulse width 
and time delay associated with one input signal, and 62 and 
©2 are the corresponding quantities associated with the other 
input signal., Thus, it can be observed that this technique 
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provides greater flexibility in setting the angular 
limits, especially the asymmetrical ones. The composite 
characteristics, however, do not enclose the fault area 
as closely as desired. 

Jackson [56] developed shaped mho and polarised mho 
characteristics, besides other shaped characteristics. 
These characteristics are shrunk along a direction 
perpendicular to the replica impedance complexor# He 
obtained these shaped characteristics with signal 
dependent phase angular limits. The signal dependency 
of the angular limits is accomplished by modifying the 
rate of fall of integrator output voltage with the appli- 
cation of an additional voltage. This additional voltage 
is a block corresponding to the anticoincidence of one 
of the normal relay inputs and an additional relay input. 
The additional relay input is proportional to the system 
current, if circular characteristics are to be shaped. 
Each additional input, along with the coincidence block 
of the normal relay inputs, is applied to a separate 
integrating circuit. Apart from the increase in the 
complexity of the circuitry, when this technique is 
employed, the shaped characteristics also do not enclose 
the fault area closely. 

Johns [53], Parthasarathy et al. [59] and Pai-Uiankar 
and Thoko [62] developed adaptive or variable pickup 



characteristics whose tripping area will be minimal under 
all balanced conditions including heavy loads and power 
swings, and automatically expand to a large area under 
unbalanced-fault conditions. Johns obtained this 
feature by adding an extra current-derived component 
to one of the relay inputs. This coitqDonent is so formed 
that it is zero under all balanced conditions and 
increases the resistive reach of the relay under 
unbalanced-fault conditions. Using the generalised 
phase comparison and the present techniques, he generated 
a variable characteristic which comprises of segments of 
restricted directional, restricted ohm and restricted mho 
characteristics. Parthasarathy et al. realised the adaptive 
feature by controlling the angular limits of phase comparison 
with negative or zero sequence voltages, arid generated an 
adaptive polarised mho characteristic, Paithankar and Thoke 
accomplished this feature by incorporating zeto sequehce 
voltage polarisation and cross polarisation in two inputs 
of a three-input double phase comparator relay, and obtained 
a variable quadrilateral characteristic suitable for earth- 
fault protection of double end infeed transmission lines. 
This technique needs the derivation of additional sequence 
components which are complex in some cases. 
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Several authors have reported about the generation 
of a quadrilateral pickup characteristic. McLaren [42] 
and Khincha et al. [48,49] have reported about its 
generation with plural arrangement of two-input amplitude 
comparators, whilst Paithankar and Ingole [57] reported 
about its generation with a multiinput amplitude comparator. 
McLaren developed a sampling technique using which any 
pickup characteristic can be obtained by comparing instanta- 
neous values of voltage and current derived at two different 
instants of time. Using this technique, he described the 
generation of parallelogram and restricted directional 
characteristics, besides some others. He explained how a 
quadrilateral characteristic could be obtained by compounding 
these two characteristics. This technique dispenses with 
the need for phase shifting and mixing circuits. 

Khincha et al. developed an amplitude comparison 
technique in which full-wave rectified, but unsmoothed, 
voltage and current combinations are compared at every 
instant during each half cycle to generate various characteri- 
stics. This technique needs no sampling circuits, delay 
circuits or amplitude-to-puls e-width converters as compared 
to the sampling technique. They described how a quadri- 
lateral characteristic could be obtained either by AND 
compounding parallelogram and directional parallelogram 
characteristics or by OR gating two directional parallelogram 
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characteristics. However, multicomparator arrangements, 
used to generate a compounded characteristic, are associated 
with time coordination problems. 

Paithankar and Ingole described the generation of a 
quadrilateral characteristic with a five-input amplitude 
comparator. The comparator is very simple. It consists of 
a switching transistor, to the base of which the rectified 
restraining input, in series opposition to the parallel 
combination of the remaining four rectified inputs, is 
applied. Trip signal is obtained at the collector of the 
transistor only when one or more of the rectified operating 
inputs exceed the rectified restraining input. Each 
operating input, in conjunction with the restraining input, 
produces one side of the quadrilateral characteristic. 
However, this scheme requires as many as five inputs and 
their rectification. Consequently, elaborate circuitry is 
needed. 

The simplest way of generating a quadrilateral chara- 
cteristic is by AND compounding restricted directional and 
restricted ohm characteristics. Humpage and Sabberwal [35], 
Hans Hoel et al.[383. Anil Kumar [50] and Ramamoorty and 
Lall [58] described the generation of a quadrilateral 
characteristic in this way. Humpage and Sabberwal employed 
two-input phase comparators, with asymmetrical angular limits 
of comparison, for this purpose. They developed a technique 
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to obtain asymmetrical angular limits. The basic principle 
involved is the determination of the coincidence of the 
relay inputs’ coincidence block and a pulse produced at a 
fixed time delay from the start of this block. The angular 
limits are a function of this time delay. The assymmetrical 
limits are obtained by providing a different time delay 
according to which input is leading the other. The demerit 
of this technique is the possibility of maloperation due to 
a spurious pulse. 

Ramamoorty and Lall developed a new two-input phase 
comparator in which assymmetrical angular limits can easily 
be derived. The comparator uses digital circuits. A novel 
feature of this comparator is the use of a digital up/ down 
counter. The counter, in conjunction with constant frequency 
pulses, counts up during the coincidence period and counts 
down during the anticoincidence period. The net count, if 
positive, at the end of one cycle time, is taken to indicate 
a trip condition. Angular limits other than + 90® are 
obtained by using pulses of different frequencies during 
coincidence and anticoincidence periods. And, the asymmetrical 
limits are obtained by providing a different ratio of coinci- 
dence to anticoincidence counting pulse frequencies according 
to which input is leading the other. This comparator is not 
affected by spurious spikes and has a uniform operating- 
of one cycle for faults anywhere in the trip region. They 
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used this comparator to produce restricted directional 
and restricted ohm characteristics, and a quadrilateral 
characteristic by AND compounding these characteristics. 

Anil Kumar [50 ] described the generation of restri- 
cted directional and restricted ohm characteristics in a 
novel way. Each of these characteristics was produced 
with a phase sequence detector fed with three appropriate 
inputs. He described two types of phase sequence detectors? 
one uses a ferrite-core logic which needs no d.c. supply, 
and the other, semiconductor flip-flops. However, use of 
either two two-input phase comparators or two phase sequence 
detectors, to generate a quadrilateral characteristic in this 
manner, is involved with time coordination problems. 

Ramamoorty and Lall [60] developed a completely novel 
way of generating a quadrilateral characteristic with only 
three input signals. The three inputs used are? IZ^, (iZj^-V) 
and (IZ^-V) where V and I are relay’s input voltage and 
current, and Zj^ is the replica impedance. One side of the 
quadrilateral is fixed, and is equal to Z^, Starting from 
the zero crossing of IZp^, three pulses, with successive 
ones adjacent to each other and of widths equal to the 
desired angles between the successive adjacent sides of 
the quadrilateral, are produced by monostable chips. 

Complete coincidence of the first pulse and (IZj^-V) as 
well as that of the third pulse and (IZ^-V) indicate a 
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tripping condition. Partial coincidence of either or 
both indicates a blocking condition. The coincidence 
of these two pairs is measured, through an AND gate by 
a digital counter, in conjunction with constant and 
high frequency counting pulses. The counter is set to 
give an output for a count corresponding to the full 
coincidence condition. 

Parthasarathy [36] developed a three-step quadri- 
lateral relay, which employed transistor circuits for 
measuring, discriminating and logic functions. It 
consists of a phase detector and a zone differentiation 
circuit. VJith four appropriate inputs, the phase detector 
generates a directional open quadrilateral characteristic. 
One of the inputs is a pulse and the detector’s operation 
is based on the extension of block-spike technique to 
multisignal phase comparison. The zone differentiation 
circuit produces three plain impedance characteristics, 
corresponding to the three zon«s of distance protection, 
either by amplitude comparison or by phase comparison of 
two appropriate inputs. The outputs of the phase detector 
and zone differentiation circuits are compounded to yield 
a three-step closed quadrilateral characteristic, 

2.3 PROPOSED RELAYING SCHEME I 

This scheme, generates a three-step quadrilateral 
characteristic. The principle of operation, hardware 
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implementation and test results of this scheme are presented 
in this section. Although the principle of operation of this 
scheme is the same as that given in the reference [36], the 
hardware implementation is different. In the present pro~ 
poSed scheme. Op Amp, monostable and digital IC chips are 
employed for measuring, discriminating, logic and other 
functions, whereas transistor circuits were used in the 
scheme described in the above reference. 

2,3.1 Principle of Operation 

•i 

This relaying scheme consists of a forward fault 
detector (FFD) and a zone detector (2D). The FFD 
produces a directional open quadrilateral characteristic 
on the complex impedance plane, whilst the ZD produces 
three plain impedance characteristics corresponding to the 
three zones of protection, as shown in Fig, 2.1, 

The FFD is basically a multi-input phase comparator 
and operates on the block-spike principle. It is supplied 
with the following inputs. 

Si « iZ^ pulse (at negative-going zero crossing of IZ^) 

$2 (IR-V) 
s V /-90^ 



Where V and I are the relay’s input voltage and current, 
is the replica impedance and R is the maximum expected 
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fault resistance. That the resultant threshold characteri- 
stics generated by the FFD is an open quadrilateral chara- 
cteristic can be explained as follows. The inputs? IZ 

Jm, 

pulse and V/ -90° , give rise to a directional characteristic 
with an operating criterion of (© -90*^) < 9 < (0 +90^) 

JLm JL 

Where 6 = arg{Z ) and 9 s= arg{v/l). The IZ_ pulse with V 

J. JL X 

and (IR-V) inputs gives rise to two blinder characteristics 
having, respectively, the following operating criteria : 
(-180°+e^f < 9 < ©> and (-180°+©^) < arg[lR-V] < 6 ^. The 

result of the above comparisons would be an open quadri- 
lateral characteristic as depicted in Fig, 2.1. 

The ZD is a set of three two-input amplitude compara- 
tors supplied with the following signals, 


jvj and jlZ^ 
jk^vj and | IZ^ 
and 1 k^vj and 


for zone 1 comparator 
for zone 2 comparator 


IZ. 


for zone 3 comparator 



where k^ and k 2 are variable settings on a potentiometer. 
The three amplitude comparators produce three impedance 


circles of radii equal to 
as shown in Fig. 2,1, 


V'l 


and 


^ a/^2 


2,3.2 Hardware Implementation 

The FFD and ZD are built with ICs and a few discrete 
components. The hardware details of the FFD are given in 
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Figs. 2.2 and 2.3. Positive blocks, corresponding to the 
positive half cycles of the signals (IR-V), V/~»90 *^ and V 
are produced at the points A,B and C respectively while 
positive blocks, corresponding to the negative half cycles 
of the same signals are produced at the points E,D and 
F respectively (vide Fig, 2,2). By means of a pulse shaping 
circuit, shown in Fig, 2,3, pulses at the negative going 
and positive going zero crossings of IZ^, are produced at 
G and H respectively. The positive coincidence of the 
blocks of (IR-V), V/ -90^ and V, corresponding to a positive 
half cycle, and 12^ pulse, produced at the negative going 
zero crossing, is measured by the combination of AND gates 
1,2, and 5, Similarly, the negative coincidence of the 
blocks of (IR-V), V/ -90^ and V, corresponding to a negative 
half cycle, and IZ^ pulse, produced at the positive going 
zero crossing, is measured by the combination of AND 
gates 3,4 and 6, The outputs of AND gates 5 and 6 are 
applied to an OR gate as shown. The output of this OR gate 
is the final output of the FFD, 

The hardware details of the 2D are given in Fig, 2,4. 
The inputs, V and IZ^, are phase split, full wave rectified 
and smoothed independently to obtain jvj'and j IZ^j respecti- 
vely. I V j is applied to a potential divider from which 
jk^vj and |k 2V | are obtained, | v | and | | are applied to 

a differential amplifier circuit to obtain | • 



IZf pulse 
f rom G 
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With a diode clipper, only positive values of (| IZ^^j “I v| ) 
are made to appear at M, Non-zero positive values of this 
quantity represent a zone 1 fault* In a similar way, only 
the positive values of{ |lZ^j - 

are obtained at N and 0 respectively. Non-zero positive 
values of these quantities indicate that the fault is within 
zone 2 or within zone 3 reaches respectively. The last two 
quantities are delayed by T 2 and seconds respectively with 
mono chips. The delayed quantities appear at P and Q. The 
outputs available at M,P and Q are applied to an OR gate. 

The output from this OR gate is the output of the ZD, which 
is available immediately for zone 1 faults and after a time . 
delay of T 2 or T^ seconds for zone 2 or zone 3 faults. 

The outputs of the FFD and ZD are applied to a final 
AND gate 7 to provide a trip signal. The block schematic 
of the entire scheme is given in Fig. 2,5. 

2,3.3 Test Results 

The relay has been fabricated and tested. Only zone 1 
static pickup characteristics have been obtained. For the 
purpose of testing, IR is taken as reference. Signals V and 
IZ_ are obtained using phase-shifting networks. First, IZ 

Jm *** 

is given the desired magnitude and phase shift. Next, V is 
given different phase shifts. For each phase shift, the 
magnitude of V is varied till the boundary of the trip 
region is identified by the operation of the relay. The . 






measuring ci rc ui ts ; F W R - Full wave rectiFier ; T 2 & T 3 Delay ci rcuits 

FlG.2-5 BLOCK SCHEMATIC DIAGRAM OF THE ENTIRE 
SCHEME I 
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relay characteristics are determined for two different 

values of at a constant value of R and also for two 

different values of R at a constant value of Z • In both 

r 

the cases. 0^ is chosen to be 80^. The results obtained 
r 

are given in Tables 2,1 and 2,2, The plots of the results, 
aloAg with the theoretical curves, are given in Figs, 2.6 
and From the results obtained, it can be observed 

that there is a close agreement between the theoretical 
and experimental curves, 

2,4 PROPOSED RELAYING SCHBffi II 

The principle of the operation of this relaying scheme 
is based on the extension of the technique reported in the 
reference [60], This relaying scheme also consists of an 
FFD and a ZD, However, the FFD of this scheme needs only 
three inputs as against four in the previous scheme. Also, 
none of its inputs is a pulse. These modifications result 
in the enhancement of reliability and also in the elimina- 
tion of the possibility of maloperation due to spurious 
spikes. In the present scheme, a pulse is used only for 
resetting a digital counter. Therefore, a spurious pulse 
will only reset the digital counter prematurely, resulting 
in the delayed operation of the relay in the event of an 
internal fault. The ZD of the present scheme generates 
three reactance characteristics corresponding to the three 
zones of distance protection. This modification leads to a 
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TABLE 2.1 

The Pickup Characteristics of Scheme I 

Q=Q0°; IZ^ = 4.0V and IZ^ =3 ,0V 
r ' r 2 


Phase Angle 
(Degrees) 

Critical value of 

V for (volts) 

^1 

Critical value of 

V for 

^2 

0 

1.5 

1.5 

10 

1.6 

1.5 

20 

1.7 

1.6 

30 

1.9 

1.8 

40 

2.3 

2.2 

50 

2.9 

2.8 

60 

4.0 

3.0 

70 

4.0 

3.0 

80 

4.0 

3.0 



The 

0 

r 

TABLE 2,2 

Pickup Characteristics o 
= 80° j IR^ = 1.5V and IR 2 

f Scheme I 

= 2.0V 

Phase 

Angle 

Critical value of 

Critical value of 

(Degr 

ees) 

V for R^( volts) 

V' for R 2 (volts) 

0 


1.5 

2.0 

10 


1.6 

2.1 

20 


1.7 

2.2 

30 


1.9 

2.5 

40 


2.3 

3.0 

50 


2.9 

3.9 

60 


4.0 

4.0 

70 


4.0 

4.0 

80 


4.0 

4.0 



— incorcticai 
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I r (volts) 
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FlG-2-6 PICKUP CHARACTERISTICS OF SCHEME I 
FOR TWO VALUES OF Zr 


Ix (volts) 



FiG.2F/ PICKUP CHARACTERISTICS OF^SCHEME I 

FOR TWO VALUES OF R 
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better fault resistance coverage in the case of faults near 
a zonal reach. The relay is built with digital ICs. This 
results in a simpler circuitry, which is another factor 
contributing to the enhancement of its reliability. 

2,4.1 Principle of Operation 

The FFD of the present scheme is supplied with the 
f ol lowing inputs • 

- IR 

$2 = (IR-V) 

$3 = V 

It yields a directional open quadrilateral characteri- 
stic as explained below. The open quadrilateral characteri- 
stic should enclose all those points, such as P shown in 
Fig, 2,8, representing forward faults on the protected 
transmission line and beyond. For P to fall v/ithin it, the 
following two conditions must be fulfilled simultaneously : 
(i) it should be above DOB and to the left of BC and (ii) it 
should be to the right of EGA, 

For the first condition to be fulfilled, the vector 
PB s= (IR-V) should be v/ithin the lagging angular range, 

0® to (ISO'^-©). This condition can be realised through 
the circuitry shown in Fig. 2,9. The monostable multi- 
vibrator, mono 1, produces a block of duration (180®-©) 
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starting from the positive going zero crossing of IR. 

Mono 2 produces a block of duration 0 starting from 
the instant when the previous block becomes zero. Only 
if (IR-V) lags IR by not more than ( 180°-©) , there will be 
complete coincidence of (IR-V) and ’©' block (the output 
of mono 2) as shown in Fig, 2,10. Hence, the output 
of AND 1 will be of full duration of 0 only if P lies 
above DOB and to the left of BC. It can be observed 
that the slope of the line BC can be changed simply by 
changing the widths of the two blocks. 

For the second condition to be fulfilled, the angle 
by which V leads IR should lie within the angular range 
of 0° to 0^, This condition is realised through the 
circuitry shown in Fig, 2,11, Mono 3 produces a block 
of duration (IBCR©^), starting at the positive going zero 
crossing of IR. There will be a full coincidence of V and 
the output of mono 3, only if V leads IR by not more than 
0^, This is illustrated in Fig, 2,l2. Hence, the output 
of AND 2 will be of full duration of (l8O°-0^), only if 
this condition is fulfilled. 

For fulfilling both the conditions simultaneously 
and thus establishing the tripping criterion, an additional 
circuitry is employed. This is shown in Fig, 2.13, The 
outputs of AND gates 1 and 2 are applied to AND 3 through 
an OR gate. The output of a clock (a standard oscillator). 
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which consists of train of pulses at a constant and high 
frequency, is also applied to AND 3, The output of AND 3 
is fed to a digital counter which is set to give a trip 
signal only if the number of pulses from AND 3 corresponds 
to the full period of ( 18(?~0^+©) , that is, only if P is 
within the open quadrilateral characteristic. Normally, © 
will be less than © • Consequently, the above measurement 
will be completed in half a cycle. Hence, with the above 
circuits duplicated for the negative coincidence measurement, 
the counter can be reset every half a cycle. Then, the 
maximum time of operation of the relay will be about one 
cycle. However, with only positive coincidence measurement, 
the time of operation of the relay will be between ^ and 1^ 
cycles. The compiete schematic diagram of the FFD is given 
in Fig. 2.13. 

The ZD is supplied with the following inputs. 


-IX^ 

and 

-IX^4-V for zone 1 
r 

-IX^ 

and 

-IX^+k.V for zone 2 
r 1 

and 



-IXj, 

and 

-IXr+k2V for zone 3 

where 

is 

the imaginary part of the replica impedance. With 


each pair of inputs, a reactance characteristic is produced 
by phase comparison prihciple. J^dr zone 1 reactance chara- 
cteristic, -IXj., (-IX^+V) and the output of a clock are applied 
to an AND gate (AND 4) as shown in Fig. 2.14. The output of 
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AND 4 is applied to a digital counter (counter 2), The 
number of pulses applied to the counter depends upon the 
duration of the coincidence of ~IX^ and (-IX^+V)* The 
level of the counter is set such that it gives a trip signal, 
only if the coincidence period of -IX^ and (-IX^+V) equals 
or exceeds 90°, This results in the generation of a 
reactance characteristic. The counter is reset by (-IX^) 
pulse. Similar circuits are used for zones 2 and 3, and the 
necessary time delays are incorporated, as shown in the 
Fig* 2,14« The outputs of the circuits of all the zones are 
fed to an OR gate. Inputs ~IX and (-IX +V), instead of 

Jm 

the normal IX^, and (IX^-V), are used to generate the rea- 
ctance characteristic. This choice makes the outputs of 
the counter 1 of the FFD and counter 2 of the ZD to be 
available at the same instant, as shown in Fig. 2.15, and 
therefore, obviates the need for special time coordination 
measures. However, special time coordination measures are 
not necessary for zones 2 and 3 only when the time delays 
are integral multiples of full cycle time. 

2,4.2 Hardware Implementation ' 

The hardware details of the FFD of this relaying 
scheme are given in Fig, 2,16, The signals IR, (IR-V) and 
V are squarewaved by means of operational amplifier (Op Amp ) 
chips. Timer 555 chips are used for building the monostables. 
The duration of the blocks, produced by the monostables, is 
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varied with the aid of (0-50 KO) potentiometers* The 
magnitudes of the outputs of the monostables and those of 
the squarewaved V and (IR-V) signals are limited to about 
+5 volts by means of zener diodes* 

The hardware details of the ZD are given in Fig. 2,17 
for zone 1. The circuitry is identical for other zones. 

The input signals -IX^ and (V-IX^) are squarewaved by 
means of Op Amp chips, and are limited to about +5 volts 
with the help of zener diodes. 

The standard oscillator is built with the timer 555 chip 
as shown in Fig, 2,18. A train of pulses, with a constant 
frequency of 15 KHz, is generated by this oscillator (clock). 
It is used for both the FFD and ZD. The digital counter is 
built with a cascaded arrangement of 4-bit binary counters 
(7493) and magnitude comparators (7485), as depicted in 
Fig* 2,19, The clear pulses required for resetting the 
digital counters are generated, just before and terminating 
at the positive going zero crossing of the concerned signal, 
by the circuit shown in Fig, 2,^, 

2*4.3 Test Results 

The relay, whose hardware details are described in 
the previous section, has boon fabricated and tested* The 
digital counter of the FFD is set to give a trip output if 
the number of input pulses exceeds 142 in a cycle. This 
corresponds to = 70® and © = 60®. The digital counter 
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of the ZD is set to give a trip signal if the no* of input 
pulses exceeds 75* This corresponds to an angle of 90°. Both 
the settings are based on the input-signal and clock fre- 
quencies of 50 Hz and 15 KHz respectively* 

Only the zone 1 pickup characteristic has been 
determined. The testing procedure adopted is as follows. 
First, the desired phase angle between V and I is obtained. 
Keeping this phase angle constant, the magnitude of V is 
decreased from a high value till the boundary of the trip 
region is identified by the operation of the relay. This 
procedure is repeated for different phase angles varying 
from 80° to 20° in steps of 5° each. Also, the nonoperation 
of the relay has been observed for a reverse fault at each 
angle, which is simulated by just reversing the direction of 
the current. A light emitting diode (LED) is used at the 
output terminals of the relay for indicating the operation 
of the relay. The results are given in Table 2,3, The 
experimentally-determined and theoretical pickup characteri- 
stics are given in Fig, 2,21, From these, it can be observed 
that there is a very close agreement between these 
characteristics. 

2,5 CONCLUSIONS 

The overview of some important distance relay characteri- 
stics, suitable for the protection of long and heavily-loaded 
EHV/UHV transmission lines presented in this chapter, leads 
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TABLE 2.3 


The Pickup Characteristics of Scheme II 
Zg = 23/0® Ohms 


Phase Angle 
(Degrees) 


IX (volts) Critical value of V 

(volts) 


20 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

70 


1.50 

4.7 

2.25 

5.3 

3.00 

6.1 

4.00 

7.1 

5.50 

8.7 

7.00 

9.8 

7.00 

9,2 

7.00 

8.6 

7.00 

8.0 

7.00 

7.8 

7.00 

7.5 

5.00 

5.3 

3.00 

3.2 


70 



Theoretical 
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O Ex pcnmental 
Zs ^ 23 [0 A 



FlG-2-21 THE PICKUP CHARACTERISTICS OF 
SCHEME II 



67 


one to the conclusion that the search for simpler techniques 
and circuitries is still necessary. Two new relays, employ- 
ing ICs and yielding a three-step quadrilateral characteri- 
stics, have been developed, fabricated and tested. The hard- 
ware implementation of both the relays has demonstrated that 
the use of ICs results in great simplicity o.f the circuits. 
Also, the test results show that there is a close agreement 
between the theoretical and experimentally-determined 
threshold characteristics, in the case of both the relaying 
schemes. 



CHAPTER 3 


FREQUENCY DOMAIN FAULT ANALYSIS OF MULTINODE 

POWER SYSTEMS 

3*1 INTRODUCTION 

Currently much interest is bestowed upon the 
development of ultra high speed protective relays. In 
order to design and develop such relays, an accurate 
determination/ simulation of post-fault waveforms or of only 
the travelling wave components, present in them, during a 
very short interval of time after fault inception, is 
essential* This necessitates an accurate modelling of 
power system components and also, the use of methods of 
computation, which are not only accurate but also amenable 
to digital computer application. The time-domain and 
frequency-domaip methods are used for the computation of 
fault currents and voltages. However, the frequency-domain 

I 

methods, though computationally cumbersome, are more ^ 
accurate [llSj. Four recent papers [3,68-70] have dealt 
with these aspects of modelling and fault analysis. 

Humpage et al* [68] have described the digital simula- 
tion of models of synchronous generators, transmission lines, 
power transformers and transducers. But, they have con- 
sidered only a simple lumped-parameter and. frequency- 
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invariant series R-L model for the transmission lines. 
Kothari et al. have presented modelling techniques for 
simulating the primary system and transducers [69] i and also 
the digital simulation of fault voltages and currents [70]. 
Although, they have considered the distributed nature of the 
line parameters, they have not incorporated the frequency 
dependence of the line parameters and the untransposed 
nature of the lines. They have derived detailed expressions 
for voltages and currents for different types of faults 
and recommended the use of bus impedance matrix for 
frequency-domain fault analysis of complex power systems. 
Johns and Aggarwal [3] have described the digital simulation 
of faulted lines with an exact frequency-domain model for 
transmission lines. That is, they have considered the 
frequency dependency and distributed nature of the line 
and earth parameters as well as the untransposed nature of 
the lines. They have modelled the line by a two-port 
[ABCD parameter] network in the frequency domain. They 
have also developed techniques for the frequency-domain 
fault analysis. However, these techniques would be appli- 
cable, with ease, only to simple systems. 

In the frequency-domain fault analysis of multinode 
power systems, the use of bus impedance matrix may not 
always be convenient in view of the fact that the bus 
impedance matrix has to be formed for a large number of 
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frequencies, and stored. With this in view, frequency- 
domain fault analysis techniques, employing bus admittance 
matrix, are developed in this chapter. Also, a frequency- 
domain Pi model of a transmission line, which is suitable 
for bus admittance matrix formation, is developed. Appli- 
cation of these techniques to determine the relaying-point 
voltages and currents for all types of faults is outlined. 
The viability of these techniques is illustrated by deter- 
mining the fault-generated components of the relaying-point 
voltages and currents for symmetrical three-phase and 
single-line-to-ground faults in tWo sample power systems. 

3,2 FREQUENCY DOMAIN. PI MODEL OF A TRANSMISSION LINE 


Under transient conditions apd in the frequency domain, 
the sending-end and receiving-end quantities, of a three- 
phase transmission line with distributed parameters, are 
related by the following equation [vide Appendix A] 
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^a,b,c 

■ : 


A B : 


V^>b, c 

R 


S3 




*^3 ^ Id 1 c 


C D 


■^»t>,c ’ 

MM. „,ri -r 


MM MM 




where 

- column vector of the transforms of the line 
to ground voltages of the phases a,b,c 
Y®,b,c _ column vector of the transforms of the line 
currents pertaining- to tfT^phases a,b, c 
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A,B,C,D - (3x3) matrices of frequency domain transmission- 
line constants 

and, subscripts S and R refer to the sending and receiving 
ends of the line respectively* 

Now, consider the Pi model of the line shown in Fig. 3*1. 
The following equations can be derived for this network • 

^a.b,c , ^ ^a.b.c ^,b,C) ^.b,c , ^a.b.o ^a,b,c 

and 

Ya,b,c _ / 7 a,b,c ^a,b,c 7 a,b,c r=a,b,c 7 a,b,Cx ■^fbjC 
ig - U2P + ^ip + ^IP ^p ^2P ^ ''r 

+ + (3.2) 

wh ere 

- three-phase transform impedance matrix of the 
series branch of the Pi model 

^IP^*^ and - three-phase transform admittance 

matrices of the shunt branches of the Pi model 
at the two ends 

and U - identity matrix of order 3* 

Comparing equations (3.1) and (3,2), and solving, we get, 

=B 

= (D - U) B~^ 

~a,b,c _ g-1 _ U) 

2P 


and 
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But, 

(D-.U)3‘'^ = (QD^ Q'“^-U)Q 

= « °m - Q S-1 

= Q \ b;i s-1 - Q b;;! s-1 (V = V 

= ® S“^, since A^ and B„ are 

diagonal matrices 

= «b;^s-1sa„ s-1-qb;;1s-i 

= A - B“^ = B“^ (a - U) 

where 

S and Q - voltage and current modal transformation 
matrices 

and -■ submatrices of modal line constants* 

Therefore, finally we have 

' = B 

and 

■^jb,c ^^^b,c _ g-1 (3^3) 

Equation (3,3) defines the parameters of the Pi model in 
terms of the A,B,C,D parameters of the line, 

3.3 FREQUENCY DOMAIN FAULT ANALYSIS USING BUS ADMITTANCE 
Mi\TRIX 

The voltage or current at any point in a power system 
network, after the occurrence of a fault, can be regarded as 
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the sum of the j^efault and fault-generated components# 

The fault analysis techniques, presented in this section, 
enable the evaluation of the fault-generated components 
in the frequency domain. The time-domain values of these 
can be obtained by any standard numerical inverse transform 
method • The prefault components can be determined in the 
usual manner by load flow studies. The sum of the pref'a'ult 
and fault-generated components gives the actual post-fault 
values. In this connection, it may be stated that the deter- 
mination of the fault-generated components separately is 
quite useful for designing and developing ultra high speed 
relays operated only with these components. As the above 
procedure invokes superposition principle, it is valid only 
for linear systems, 

3,3,1 Application of the Superposition Principle [3] 

The prefault three-phase power system transform network 
is shown in Fig, 3, 2(a), where only the fault node is shown in 
detail. The open three-pole switch represents the unfaulted 
condition. It can be replaced by a three-phase transform 
voltage source » having polarities as shown in Fig, 

3, 2(b), are the transforms of the three prefault 

phase voltages at the fault point. A fault can be simulated 
by connecting another three-phase transform voltage source 
Vpp^’'^ in series opposition to as shown in Fig, 3,3, 

Choice of will depend upon the type of fault to be 




FIG-3-2 PREFAULT NETWORKS 



FIG.3-4 NETWORK FOR EVALUATING ONLY THE 
FAULT GENERATED COMPONENTS 
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simulated. For example, for a symmetrical 

three-phase fault, and for one-phase-to-ground fault 
say on phase 'a’ , = V^p* The values of ^p and*^P 

in this case (i.e. one-phase-to-ground fault) are such 
that the fault currents in phases b and c are zero* 

Whatever the type of fault, the fault resistance network 
remains the same. The presence of unwanted resistances 
is rendered ineffective by setting the currents in them 
equal to zero. Application of Superposition theorem to the 
network of Fig. 3.3 yields the network of Fig. 3.4. The 
latter network is used for determining the fault-generated 
components. 

3.3.2 Representation of the Fault Node 

A fictitious bus, called fault bus hereafterwards, 
is created at the point of fault. At this bus, a series 
branch of the three-phase voltage source and the 

three-phase fault resistance Rp * exists to ground, as 
explained previously, where, , 

na,b,c 

rtp 

For the convenience of nodal analysis, it is replaced by 
an equivalent circuit consisting of a 3-phase current 
source and a 3-phase conductance Gp*^*^ in parallel, as 


R_ 


R 


g 


R. 


R4r+R-„ 

f g 
R„ 




(3.4) 
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shown in Fig, 3.5. For equivalence, it can be readily seen 
that. 


"rdjbjC 

■Jp 


(3a,b.c^.b,c 


.a,b,c rna,b,c-i-l 


(3.5) 


For salc« o-f gfenerali^^ Xet 


Z a,b,c 
F 


Qa,b,c 

Rp 


X7a,b,c 


pa,b,c 

^F 


(3,6) 


3.3.3 Nodal Equations and Bus Admittance Matrix 

Assuming (N+l) number of buses, including the fault bus 
F, and with ground as the reference node, the followir^ 
frequency-domain nodal matrix equation can be developed for 
a power system network to determine the fault-generated 
components (vide Fig, 3,4). 

V® » ^ j ® ^ ^ ^ c I V® » h , c 7»a , b , c n 

\2 ••• ^IN , ^IF ^ 

^a,b,c ^^,b,c ^na,b,c i ':^,b,c ?^,b,c.. q 

^21 ^2 . ••• j ^2F ^2 - ^ 

• • • 1 * * ~ * 


«a,b,c «|,b,c , ^a,b,c ^,b,c '^,b,c 

^N1 ^2 ••• ^NN \ W ''N 
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where 

» - column vector of the transforms of the 

voltages of phases a,b,c at ith bus for all 

•••> NjF 



It can easily be recognised that ^^ 3 '^ is the bus 
admittance matrix. Equation (3,7) or (3.8) can be solved 
for determining any desired bus voltage transforms. Knowing 
the bus voltage transforms, the current transform in any 
element can be evaluated. However, in order to solve the 
above equations, should be known. Since is 

related to through eqns, (3,5), it is necessary to 

determine Vp^^*^ for different types of faults. 

3.3.4 Determination of for Different Types of Faults 

(i) 3-Phase Symmetrical Fault 

An examination of the post- fault power system network, 
shown in Fig, 3.3, reveals that should be equal to 

^a,b,c, whi^ch is tantamount to closing the 3-pole switch 
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of Fig. 3,2(a), for simulating this type of fault. Hence, 
for this type of fault. 


^a,b,c ^ "^jbjC 
VpF - Vpp 


(3.10) 


(ii) One-Phase-to-Ground Fault on Phase 'a’ 

This fault is simulated by letting = \^p, and by 
letting V^p and Vpp to assume such values that the fault 
current transforms in phases b and c are zero i.e,, 

Tp = 1[p = 0. The values that ^p and \^p should have, in 
order that Ip = Ip =0, can be determined as outlined 
below. 

Solving equation (3.8), we get 


r«a,b,c/^,b,c^-l ^,b,c __ •^,b,C’j^a,b,c _ ^a,b,c 
LYp ) % ^FF J'^F - '^F 


r^,b,c ^,b,c 
F FF 


(3.11) 


Also, at the fault bus, the three-phase fault current 
transform is given by (vide Fig, 3,5) 


ia»bfC _ -jatbjc ^ "^,b,c Ya,b,c 


TrafbjC ^,b,c , Ya>b,c T^,b,c 

“ ^FF ^ ^F ^F 

Solution of equations (3.11) and (3,12) gives 


(3.12) 


FF 


^a,b,C(:^,b,C)-l :^a,b,cj-l^-ja,b,c 


n 


i.e.. 


(3.13) 
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where 


^a,b,c ^ ^,b,c _ [^,b,c _ :^,b,c ^ ^,b,C(^.b,C)-l ^.b.cj-l 

(3.14) 


and 




Now, let 


^aa 

^ab 

^ac 

2ba 

^bb 

2bc 

jca 

l: 

^cb 

fCb 


^a,b,c 


Solution of equation (3.13), with the condition 


(3.15) 


?b TC 


•p “ -^p 


I„ = 0, results in 


^,b,c 

VpF 


^aa 


jba 


gca 




(3.16) 


(iii) Phase- to-Phase Fault between Phases b and c 


For the simulation of this fault, should be suchi 

that the following conditions are satisfied. 

^F "" ^FF ~ ^PF *" ^PF 

T® = 0 (3.17) 

+ Ip =0 


and 
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Solution of equations (3,13) and (3.17) yields, 

0 
1 

-1 




- v£.„ 

PF ^PF 

(2bb^^bC_^Cb^^CC) 


(3.18) 


(iv) Double-Phase- to-Ground Fault Involving Phases b and c 

This fault is simulated by letting = V^p, V^P = "^p 
and by letting \^p to take up such a value that Tp = 0. 

Substitution of these in equation (3.13) leads to the 
solution. 


vIp = [z®'’ z®‘=] 


'FF 


2bb 

^bc 


VpF 

zcb 

^cc 


^PF 


(3.19) 


, Bolted faults are simulated by assuming a very small 

/ 

value for Rp and by taking = 0. 

3.4 INVERSION TO THE TIME DOmiN 


After having determined the frequency-domain values of 
tho relaying -point voltages and currents, it is necessary 
to determine their time-domain values. For this purpose, 
the modified Fourier transform techniques, described in [ll6], 
are employed. The inverse modified Fourier transform, given 
in equation (3.20) below, forms the basis of the method by 
which the wide frequency spectrum of the disturbance is used 
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to determine the corresponding time variation of any 
voltage or current of interest* 

f(t) =S2||ail. J 7^ (a+J «) exp(j Mt)d« (3.20) 

where 

» - angular frequency as well as the transform 
parameter 

a - frequency shift constant which is used to reduce 
aliasing errors when the integral is evaluated digitally 

and "f^^Ca+j ») - modified Fourier transform of the function, 

f(t). 

For the numerical evaluation of the above integral, it 
is necessary to truncate the infinite range of the integral 
to some finite value, say (r n» This truncation gives 

rise to Gibb's oscillations, vhich are suppressed by using 
the Lanczos Sigma factor, defined by 

^ ° (3.21) 

With these modifications incorporated, and after some 

manipulation [116], equation (3*20) becomes 
■ , . Q 

j gp cosMt dw (3.22) 

o 

where P = Re[fj^(a+j «»)] 
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Normally, the value of the function is required at 
several instants of time. In such a case, s can be 
computed beforehand and stored. 

3.5 EVALUATION OF THE TRANSMISSION LINE PARAMETERS 

The capacitances and shunt conductance are voltage 
dependent, whilst the inductances and series resistance are 
frequency dependent. The general effect of corona discharge 
is to increase the capacitances and the conductance to earth. 
There is no quantitative information available at present 
to determine this effect in the case of multiconductor 
lines [ll6]. Accordingly, the capacitances are assumed to 
be independent of voltage in the present work. The 
conductance to earth is neglected. The self and mutual 
capacitances are determined by the potential coefficient 
method. 

The frequency dependent series self and mutual 
impedances of the phase conductors and earth wires are 
determined by using the formulae given in reference [ll6]. 
Although these formulae permit the consideration of a 
stratified earth, the earth is considered only as a single 
layer in this work. The presence of earth wires in the 
final mathematical formulation is eliminated by math^atical 
manipulation with the assumption that the potential of the 
earth wires is zero. A computer programme is developed for 
the above. It is used as a subroutine in the main programme 
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which is used for the computation of relaying-point 
voltages and currents. 

3.6 NUMERICAL EXAMPLES 

The application of the fault analysis techniques 
developed is demonstrated by considering two sample power 
systems, and by determining the relaying-point voltages and 
currents for 3-phase and one-phase-to-ground faults. 

3.6wl Three-Bus Power System 

The 3-bus power system, given in reference [117] and 
reproduced in Fig. 3.6, is considered with the following 
change. A 400 KV quad conductor three-phase single-circuit 
transmisssion line, whose tower configuration and conductor 
spacings are given in Fig. 3.7, is considered instead of 
that given in the reference. An earth resistivity of 
100 ohm-m is assumed. Other basic data of the line are 
given in Appendix B. Identical (generator) sources are 
considered at buses 1 and 2. The data of these sources 
are given in Appendix C. The faults are simulated at bus 3* 
Under these conditions, the one-line diagram, of the equivalent 
network of the system, suitable for computing the fault- 
generated components, would be as shown in Fig. 3.8 where 
^I’^l^^l’^l matrices of the frequency-domain trans- 

mission line constants for a length of 200 km and A 2 to D 2 » 
for a length of 100 Km. The bus admittance matrix can be 
written by inspection as shown below. 




^«-6'78m-H 


^*4 9 K“ 
0-3048m 


FIG.3-8 ONE-LINE DIAGRAM OF THE EQUIVAL 
NETWORK OF 3-BUS POWER SYSTEM 
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-1 
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-B, 


-1 


-b: 


2B2^A2+Yp*^^*‘^ 


(3.27) 

The frequency-domain voltages and currents at the relaying 
point, considered at bus 1 on the line connecting buses 1 and 
are determined froift five'- following equations, which are obtained 
from the solution of the nodal equations of the system. 


va»b,c _ [2B2^-(2B2^/^+^*^’^)B2(B^^A^+B2^A2 +'^'^'^-b5[^)3"^ x 


v^*b,0 ^^,b,c 

T TF 


(3.28) 


and 




(3.29) 


’^fh,c determined, as outlined below, in order to compute 


^a,b,c f^rora eqns. (3.28) and (3.29) respectively, 
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(i) 3-phase symmetrical fault. 

For this fault, 

[vide,eqn. (3.10)] 


V^[s sin9 ^+ m^cos(Pq3/(s^+ w^) 

Vjs sin( 9^-120°)+ cos( 9^-120° ) 3/(s^+ w^) 

Vjjj[s sin(9Q+l20°)+ cos ( 9^+120®) ]/(s^+ 


(3.30) 


where 

Vm - peak value of the prefault voltage at the fault bus, 
- nominal anguldr frequency of the system, 

9 q - fault initiation angle 

and s .'-i-' (a+3 w ), the modified Fourier transform parameter. 


(ii) One-phase-to-ground fault on phase 'a*. 

For the sample system, the matrix as defined 

by eqn.. (3,14), is given by 

v^a,b,e ^ ^,b,c ^ [2B2^A2-2B2^(B*^A^+B2^A2 +y|'^'‘^-B“^)bJ^3“^ 

(3.31) 

Using eqn. (3.16) in conjunction with the above matrix, 
is determined. 
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Assuming as equal to 1.0 p.u. for both the types 

of faults and employing inverse transform techniqueSf des- 
cribed in Section 3.4, the time— domain values of the relaying- 
point voltages and currents are evaluated. For this purpose, 
the following values, for the parameters associated with the 
inverse transfoirmation, are chosen. 

n =» 60,000, a - 750, n =: 500 and time step = 0.05 ms 

The results obtained, for = 0° and = 100 Ohms, 
are plotted in Figs. 3,9 and 3*10, with the phase voltage and 
line current bases chosen as equal to 400V2/V'3 KV and 
600 KA respectively, 

3.6,2 Four-Bus Power System 

The 4-bus power system, given in reference [3] and 
reproduced in Fig, 3,11, is considered without any modifica- 
tion, The tower configuration, conductor spacings and other 
basic data of the transmission lines of the present system 
are the same as those of the previous system and also as 
given in the above reference. All the lines are of 128 , Km 
length, and an earth resistivity of 100 Ohm-m is assumed. 

In this case, the sources are represented by fault levels at 
the corresponding buses [3], The fault levels are converted 
to equivalent circuits as explained in Appendix D. A ratio 
of zero sequence impedance to positive sequence impedance, 
equal to 0,5 at 50 Hz, is assumed for this purpose. 
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The faults are assumed to occur, at a distance of 

25.6 Km from bus 1, on the line connected between buses 1 
and 2, Under these conditions, the one-line diagram of 

the equivalent network of the system, suitable for computing 
the fault generated components, would be as given in Fig. 
3.12, where to D^, A 2 to and A to D are the matrices 
of the frequency domain line constants for lengths of 

25.6 Km, 102.4 Km and 128.0 Km respectively. The bus 
admittance matrix can be written by inspection as shown, 
in eqn, (3.32), given on the next page. 

Solution of the nodal equations of the system gives the 
voltage and current transforms at the relaying point which 

A 

is considered at bus 1 on the. line connected between buses 
1 and 2 . They are given below, 

(3.33) 

where 

M, = (A + B = 

1 

(A + 0,5B 

and M 3 « (SB-^A + B'^A^ + 

And, - B-^ (A^VJ»^»" - 


(3*34) 
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where s ^ '«a*b,c 

^ 4' 1 1 *p Vpp j 

*'TrQ, Id ' c . . 

* » required in the equations above, is determined 


as outlined below. 


i) 3-phase Symmetrical fault 


Vp|^»^ is given by eqn. (3,30) 

ii) One— phase- to-ground fault on phase ' a’ • 
The matrix ig given by 




(3.35) 


i-l 


where 

Mg a [B3_(M3-B’’%^-2B''^)B^}' 

As in the previous case, is computed using 

— 1 . ^ 

eqn* (3,16) in conjunction with the above matrix, Z®' * , 

For both the faults, V is assumed to be 1*0 p.u,, 

luaX 

Rj as 0,00001 and 9 ^ = 90°, The time-domain values are 
obtained by using the following values for the parameters 
associated with the inverse modified Fourier transformation^ 




* 10^, a a 1000*0, n a 500 and At a 0*05 ms- 

The results ^’’obtained' 

;l ' ' i jib'll 


3*14, with the bases chosen as 


. iij,a 


, 3,13 and 


^i'betow 


. ' r , 




KV 


Line ;curreht ' 



Phase vcltcges ana =.nc currents 'n 



Time in ms 


‘■J);?.'-.' I’i'i a " 

f»0NENT5 


F 16-3^13' 





Time in ms 
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3 •6. 3 Computer Programme 

Based on the procedure outlined in the earlier 
sections, a computer programme is developed to solve for 
the time domain values of the relaying-point voltages and 
currents* A simplified flow chart of the programme is 
given in Fig* 3*15* The steps involved are briefly out- 
lined below, 

1* Read the data of the transmission linesy sources, fault 
resistances, prefault voltages at the fault bus, fault 
Initiation angle and inverse modified Fourier transfor- 
mation parameters, 

2, Calculate the self and mutual capacitances of the 
transmission linos through a subroutine* 

3* Calculate the P^s (k » 0,1,2,,,,, n-l)* as defined in 
eqns, (3.24) and (3,25), for voltage and current 
transforms separately, using the following steps. 



i) Initialization, i,e,, set w =0*0* 

ii) Set the transform parameter, s = (a+j«), ' 

iii) For » = 0,0, assign w ss 10,0 , 

iv) Calculate the transmission lino's frequency-dependent 
self and mutual sorios impedances through a subroutine# 

v) Calculate the resistances (R)> land Inductances (L) of the 
line from, thO:a'bo,ve;i.;::,':,:.;MU‘;r.in 


When 
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vii) Calculate the per unit length series impedance 

and shunt admittance of the line as 

a function of $. 

viii) Calculate the eigenvalues and eigenvectors of the 

matrix through a subroutiner 

ix) Calculate the voltage and current modal transformation 
matrices i S and Q* 

x) Calculate the frequency-domain A, B,C,D parameters of 
the line, 

xi) Calculate the source admittances as a function of s ► 

xii) Calculate the relaying-point voltage and current 
transforms using appropriate equations. 

xiii) Calculate P^s for the above# using equations (3r24) 
and (3*25) * 

xiv) Increment ta by Au « Q /n and go to step 3 (ii)* 

4* Using eqn# (3*26)# compute the time-domain values of 

the voltages and currents at all the desired instants of 
time. ,, 

The steps 3 (iii) and 3 (vi) are necessitated by the 
fact that the frequency-dependent series self and mutual 
impedances of the line cannot bo evaluated at » = 0.0 by 
the formulae given in reference [1163<* Hence, they are 
evaluated ati an angular frequency cjlp,sq to zero* Subroutines 
have boon employed also for complex matrix inversion as well 
as multiplication* 
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3*7 APPLICATION OF FAST FOURIER TRANSFORM 

The fast Fourier transform (FFT) technique, which 
requires much less computation time, was first developed 
by Cooley and Tukey [l2l3 and later, its applications were 
discussed by Cooley, Lewis and Welch [122]* Application 
of this technique to the computation of electrical transients 
was described by Ametani [l23], who showed that considerable 
saving in computation time could be achieved in this way. 
However, it is shown that the application of FFT to the 
computation of power system transients, especially those 
under short-circuit conditions used for assessing or 
predicting the performance of high speed and ultra high 
speed protective relays, does not ^ eld as much computational 
benefit as was reported [123] • It is also shown that the 
application of this technique imposes certain restrictions* 
Accordingly, this technique is not used in this work while 
the time-domain values are computed. 

3*7.1 Problem Formulation 

In order to apply the FFT algorithm, equation (3,26) 
has to be rewritten as given below, 

n-1 ; 

f(r. At) = exp(a* A t,r) Re [ exp(j ~ kr)]| 

k=0, 

r * 0, 1,2 , ... , (m— l) 


(3.36) 
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where 


n = number of frequency-domain samples 


and m =* number of time-domain samples. 


The above equation is modified as follows : 


fix* At) 


Sf“e.p(a. .t.r).Re[2:;i^ 

i==0 {^k=0 

exp(j 2S fcr)! X exp(j 2$ n)] 


*= ^“exp(a. At.r),Re[^ B^^^,exp(j ^ ir)3f 

r — Of If **ft (j®*"l)» (3«37) 

where 

. m-1 

bU) exp(j % kr) » r = 0,1, (m-l) 

(3.38) 

and p sa n/m, an integer. r | * 

First, are calculated by using the prograimne 

* -T 

given by Ametani [1233* Next# f{rt At) is calculated for 
r = 0,1, (m-l) by using equation (3.37). The following 

observations can be made with regard to the application of 
the FFT algorithm. 
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1* Complex arithmetic is involved because of the terms» 
exp(j ^ kr) and exp(j 2a ir)* 

2, For the most efficient utilization of the algorithm, m 
should be equal to 2^ wher# q is an integer. This 
imposes a restriction on the number of time steps, 
which can be chosen, if the algorithm were to be 
utilized most efficiently, 

3. There is a blanket restriction on the choice of the time 
step At, as explained below. 

Comparing equations (3. 26) and (3.36), we get 

4. . or (3.39) 

The choice of Q depends on the problem on hand [ll6]. 
Therefore, once a value is chosen for q , then A t is 
fixed by virtue of equation (3.39). If we are interested in 
computing f(t) at intervals of (2, At) only, then we \^ould 
be unnecessarily calculating double the number of values o^ 
f(t). On the other hand, if we use equation (3.26) directly 
for the computation of f(t), then we can choose any conve- 
nient and necessary value of At. 

3.7.2 Numerical Example and Results 

For the sake of comparison, a numerical example is 
worked out, with the aid of the digital computer DEC 
System 1090 at I.I.T., Kanpur* to determine separately the 
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"time required for the computation of the frequency- 
domain values and that required for the computation of 
time-domain values when equation (3,26) is used (direct 
method), and when equations (3.37) and (3,38) are used 

4 

(FFT method). For this purpose, a 200 Km, 3-phase, single- 
circuit, 400 KV, quad-conductor transmission line, inter- 
connecting two large power networks, is considered. The 
tower configuration, conductor spacings and other basic 
data of the line are given in Fig, 3.7 and Appendix B* 

The earth resistivity is assumed as equal to 100 ohm-m. 

The sending and receiving end source networks are assumed 
as resistive, having resistances equal to l/9th and l/4th 
respectively of the real part of the appropriate modal 
surge impedances. The frequency dependency of the line 
parameters and the un transposed nature of the line are not 
disregarded. For a singi e-line- to-ground fault on phase 'a’ 
at the receiving end of the line with the assumption that the 
prefault peak voltage at the fault point is 1,0 p,u,, the 
frequency-domain values of the fault-generated components 
of all the phase voltages and line currents at the sending 
end of the line for various desired frequencies are computed, 
multiplied with the sigma factor and stored. Then, the time- 
domain values are computed for various conditions, using the 
following data* 

n = 512,1024.1 n « 62831,8 rad/s> a = 750 

m = 32,64,128,256,512,1024 
and number of functions evaluated =1,3,6, 
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The computation time taken has been determined by 
running the programmes on the DEC System 1090 computer 
at Kanpur# and the results are given in Tables 

3*1 and 3,2* 

3,7.3 Discussion on the Results 

The results given in Table 3,1 show that the FFT 
algorithm is not as fast as it has been claimed. For 
example, when n = m = 1024, the theoretical value of 
the ratio of the computation times needed by the direct 
method and tho FFT method would be 51 as given by Ametani, 
whereas it has come to be only 21,507 in this work. The 
reason for this difference is that Ametani has employed 
the complex equivalent of equation (3,26) in his studies. 
However, the equation (3,26), which involves only real 
arithmetic, can easily be rewritten in a form, suitable 
for the application of FFT, as described in Section 3.7,1. 

From the results given in Table 3,2, it can be 
observed that the relative computational relief provided by 
the FFT algorithm would be only marginal when the number of 
time-domain samples required is small. Generally, either 
for designing or for predicting the performance of high speed 
and ultra high speed protective relays, the time— domain 
samples are not needed for more than 1 cycle. And, with a 
CVT of cut-off frequency of 2,5 KHz, the sampling interval 
should be at least 200 ps. Therefore, the sampling interval 



TABLE 3.1 

Comparison of computation time in seconds taken by (i) direct method (t^) and (ii) 
using the FFT algorithm (t 2 ). of frequency samples = 1024 
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TABLE 3.2 

Comparison of the computation times required for the evaluation of the frequency- 
domain and time-domain samples 

Number of functions evaluated 6. 
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of 100 p.s» chosen in the present numerical example^ is 
more than adequate. With this sampling rate, the number 
of time~domain samples (ra) required for one cycle would be 
200. Under such circumstances, the FFT algorithm provides 
only marginal relief. For example, for n = 512 and m = 256, 
the time required for the computation of frequency— domain 
samples was 341.929 seconds, whilst that required for the 
computation of time-domain samples was 16.916 seconds by 
the direct method, and 4.352 seconds by the FFT method. Thus, 
the overall computation time has decreased from 358.845 to 
346.281 seconds only. 

3.8 CONCLUSIONS 

In this chapter, two new contributions have been made? 
the frequency-domain fault analysis of rault|.node power systems 
using bus admittance roatrixi and the development of a frequency- 
domain Pi model for a transmission line, suitable for this 
analysis. The viability of the fault analysis techniques has 
been demonstrated by employing them to the computation of 
relaying-point voltages and currents in two sample power 
systems. Although the sample power systems considered are 
relatively simple, the application of these techniques to 
larger systems would be a matter of straightforward extension. 

In the nodal method using bus impedance matrix, the bus 
impedance matrix has to be formed for a large number of . 
frequencies and stored. In cases, where this imposes 
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restrictions on computer memory, the analysis developed 
in this chapter can be employed. The computation of 
relay voltages and currents is hardly needed over more 
than half a cycle in developing ultra high speed relays. 
Consequently, the frequency-domain method of computation, 
which is most accurate, can be employed without undue 

computational effort being necessary. An accurate model 

i 

for the transmission lines has been used. And, the 
choice of source models has been based on the practices 
available in the literature. 

The application of the FFT technique to the computation 
of transients, in electrical power systems under short circuit 
conditions, does not provide, appreciable computational relief 
when the fault data is needed only for a short period. The 
results of a numerical example have confirmed this fact. 



CHAPTER 4 


TRAVELLING WAVE RELAYING OF EHV/UHV TRANSMISSION 

LINES 


4.1 INTRODUCTION 

The ultra high speed (u.h.s.) clearing of faults improves 
the transient stability of a power system because the tran- 
sient rotational kinetic energy introduced into the system 
is proportional to the square of the fault clearing time 
[72], As a matter of fact, the u,h,s. clearing of faults at 
any location, particularly on ehv/uhv Intertie lines and 
close to large sources of generation, will reduce the system 
acceleration more than any other form of dynamic control 
which can be used only after the system is being accelerated. 
The fault clearing time depends on the speed of both the 
protective relay and the associated circuit breaker. With 
the development of u.h.s. circuit breakers in the recent 
past [1,2], attention is now focussed on the development 
of u.h.s. protective relays. The development of u.h.s. 
protective relays has been possible with the utilization of 
the travelling wave phenomena. 

Quite a few authors have reported about the development 
of u.h.s. protective relays based on travelling wave pheno- 
mena. After the occurrence of a fault, the voltage and 
current at any point in a power system can be regarded as the 
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sum of pi'efault and fault-generated components# The fault- 
generated components manifest themselves as travelling 
waves during the first few milliseconds after the occurrence 
of a fault* The relaying schemes developed so far, utilize 
these travelling waves for fault detection, 

Takagi et al, [80, 8l] and Akimoto et al* [82,83] 
defined a relaying signal, which is formed from the sending 
and receiving end voltages and currents and also, the surge 
impedance of the line. This relaying signal is zero if no 
internal fault exists* and is finite if an internal fault 
exists. For the protection of three-phase lines, the modal 
components of voltages and currents, and appropriate modal 
surge impedances are used. Also, three relaying signals, 
one for each mode, are employed. This is a differential 
protection scheme. Its demerit is that it requires the 
exchange of quantitative information between the ends of 
the protected line. 

They have also proposed two fault location schemes. In 
one scheme, a discriminant, which is a function of sending 
end quantities at two different instants of time and also, 
the surge impedance of the line, is used. In the other, the 
discriminant is a function of both the sending and receiving 
end quantities. In both the schemes, the fault location 
involves the determination of the two time instants which 
would make the discriminant zero. The time interval between 
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these two instants is a function of travel time of the 
waves between the relaying and fault points. But, these 
schemes are too involved to be easily realisable. 

Yee and Esztergalyos [72], Chamia and Liberman [73], 
Esztergalyos et al. [74], Matele [75] and Carter [76] have 
reported about the development, design, application and 
operational experience of an u.h.s. relaying scheme. It 
has been developed by ASEA, Sweden, and tested on BPA's 
(Bonneville Power Administration, U.S.A.) 5CX) KV system. 

It needs the exchange of only qualitative information between 
the ends of the protected line. Consequently, it places 
loss demanding requirements on the communication channel. 

It is basically a directional comparison scheme. In this 
scheme, the fault-generated components of the voltage and 
current of the first travelling waves reaching the relaying 
point are used to determine the direction to the fault. 

These components are of opposite sign for a fault ahead of 
the relaying point, and are of the same sign for a fault 
behind. Thus, only for an internal fault, they are of 
opposite sign at both the ends of the protected line. 

Vitins [71] has developed a method, by means of which 
the location of the fault can be determined through a time 
delay equal to twice the travel time of the travelling waves 
between fault and relaying points. He has elaborated the 
application of a correlation technique for the determination 
of this time delay. 
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Johns [77] and Johns and Aggarwal [78] have reported 
about the development of an u«h«s« directional comparison 
relaying scheme* They have employed three pairs of relaying 
signals, one for each propagation mode* The signals of each 
pair are formed from the corresponding modal components 
of the fault-generated components of relaying-point 
voltages and currents. The direction to the fault is 
determined by the sequence in which the relaying signals 
of each mode exceed a preset level. The qualitative infor- 
mation about the direction to the fault is exchanged through 
a carrier communication channel. 

Vitins [85] has proposed a similar approach for deter- 
mining the direction to the fault. In this scheme also, 
the direction to the fault is determined by the sequence 
in which two relaying signals exceed a threshold level. 
However, there are two basic differences. A geometric 
approach has been employed to determine the above sequence. 
The three-phase protection has been provided with three 
signal pairs as in the previous scheme. But, the signals 
of a pair are formed from the fault-generated components 
of a phase voltage and the corresponding line current, in 
the present case. 

Ramamoorty and Verma [84] described a relaying principle 
in which a discriminant is used to determine the direction 
to the fault. The discriminant is a function of fault- 
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generated backward-wave voltage and current* It exceeds 
a preset value for an internal fault. 

With a view to provide simpler and reliable alternatives, 
two new relaying schemes, which utilize the characteristics 
of the travelling waves for detecting faults, are developed. 
One scheme is based upon the amplitude comparison of two 
relaying signals* The other scheme utilizes the reflection 
properties of the travelling waves to provide the fault 
locating capability. The relaying principles, involved in 
these two schemes, are validated by the results obtained 
from the digital computer simulation of faulted sample e,h*v* 
power system networks* 

4.2 AMPLITUDE CX)MPARISON FffiLAY SCHEME 

In this section, a relaying scheme, based upon travelling 
wave phenomena, whose operation depends upon the amplitude 
comparison of two input signals, is developed. The results 
of the digital simulation of faults in a 4-bus system, which 
confirm the validity of the relaying principles involved, 
are presented* 

4*2,1 Principles of Operation 

The underlying principles of the relay can be best 
explained by reference to a single— phase transmission line 
interconnecting two large power systems as shown in Fig, 4*1, 
Their extension to three-phase lines is explained in the next 
section. After the occurrence of a fault, the voltage and 




G-4'1 A TRANSMISSION LINE INTERCONNECTING 
TWO POWER SYSTEMS 




FIG-4-2 LATTIG 
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current at any point in the system can be regarded as the 
sum of prefault and fault-generated components, that is. 


V = Vpf + Vf 

and 

i = ipj + i^ 


(4.1) 


where the subscripts pf and f refer to the prefault and 
fault— generated components respectively. 

For a fault on a transmission line, the fault- 
generated components manifest themselves as travelling waves 
and travel towards both the ends of the line. Assuming a 
lossless line, for simplicity, the fault-generated components 
of voltage and current, at any point on the line, are given 
by the solution of the well known wave equations as, 

v^ * fj^(t-x/u) + f2(t-fx/u) 

if - ^ [f.(t-x/u) - f^Ct+x/u)] (4.2) 

where 

f^ and f2 - arbitrary functions determined by the 
boundary conditions 

t - time from the instant of fault occurrence 
u - velocity of propagation of travelling waves 
- surge impedance of the line. 


I 



From the above equations (i.e. eqn. (4.2)), it can 
be seen that f^(t-x/u) represents a forward travelling wave 
aiiu f 2 (t+x/u), a backward travelling wave. Also, the 
forward travelling waves of voltage and current are of the 
same sign, whilst the backward travelling waves of voltage 
and current are of opposite sign. 

Thus, the travelling waves of the fault-generated 
components of voltage and current, originating at the fault 
point F, travel towards the ends of the line where they 
undergo reflection and refraction. If the voltage refle- 
ction coefficients at the sending and receiving ends are 
p and p_ respectively, these events can be depicted by 

OX 

means of a Lattice diagram as given in Fig. 4.2, where '^ff 

is the fault generated component of the voltage at the 

fault point. For example, if the fault is a bolted one, 

v^^ will be equal to the prefault voltage at the point 
f f 

of fault. 

The proposed scheme makes use of a signal pair of the 
form given below for each of the relays at the ends of 
the protected line. For the sending end relay, the input 
signals (in p.u.) would be 
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where 



fs 


and 


fault-generated components of voltage and 
current in p.u, at the sending end 


^set *" impedance of the line in p,u. 


For a fault on the line, 
|S 2 | at various instants. 


the values of | | 

deduced from the Lattice diagrafti 


and eqns, (4.3), are given in Table 4.1, is the tfavel 


time of the waves to travel from the fault point to the 


sending end. Since 1 Pgl fa always less than 1,0, except 


when the sending end is on open circuit, j S^j >|S 2 jfrom 
till another wave pair arrives there. 


Next, consider a fault on any of the lines feeding 
the sending end bus. One pair of the travelling waves 
generated on this line travels towards the sending e^M bus 
and undergoes both reflection and refraction at the bus. Due 
to the refracted waves, we have 


Pt Vfjf and 




where is the refraction coefficient and Z^, the surge 
impedance of the faulted line. Then, the relay signals are 
given by 
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Si - “t ''ff - Kset "t ''ff 


^ff 

S2 = Pt Vfj + 

0 


Ot ''{f 


(4.4) 


Therefore, it can be seen that < 32^' 


till 


these waves return after reflection at the remote end of 
the protected lino. Summarising, we have for some duration, 

S^l ^ I ^2 1 for a fault ahead of the relaying point 
and 

I S^j < for a fault behind the relaying point (4,5) 

An identical relay arrangement is provided at the 
receiving end. If ■ ^ ^2| both the ends, an 

internal fault is indicated. The information about this is 
exchanged between the ends of the protected line through a 
fast communication channel. If one end of the' line is on i 

open circyit, then the relay at the other end iis arranged 
to act independently of the communication channel information. 

4.2,2 Extension of the Principles to the Protection of 
3-Phase Lines 


A distributed-parameter 3-phase transmission line can be 
described in the frequency domain by a set of decoupled modal 
differential equations given below [vide Appendix A] 

.2«(m) 

2rT(m) 

. --TS Ym ^ 


,2*(m) 


r? m = 1,2*3 


(4.6) 
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where 


Ym 


- mth mode’s components of phase voltage 
and line current transforms 

- propagation constant of mth mode. 


These eguations show that the wave propagation in a 3— phase 
line can be considered in terms of three separate and 
independent components, each possessing its own propagation 
constant and surge impedance. The three separate and 
independent components are called modal components, and are 
obtained by modal transformation, using the voltage and 
current modal transformation matrices S and Q, defined in 
Appendix A. 


It has been the practice to set the protective relays 
with the assumption of an ideally transposed line. With 
this assumption, S and Q become equal and frequency inde- 
pendent. In addition, any convenient value can be chosen 
for them, provided that they give rise to decoupled 
differential equations. In this work, the Karrenbauer modal 
transformation matrix [llS], defined below, is used. 


1 1 1 
1-2 1 
11-2 


(4.7) 


Also, with the same assumption, the modal surge impedances 
of the line are given by. 
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and 

^0^^ " ~ modes 2 and 3 

(4.8) 

whore 

mutual series impedance transforms of 

the line 

Y_#Y_ - self and mutual shunt admittance transforms of 

a ill 

the line. 


The results of Section 4,2.1 apply equally well to 
each of the three modes of propagation of travelling waves 
described above. Accordingly, three pairs of signals, one 
for each mode, are employed. For the sending end relay, 
these are given by 


c(m) _ Jm) p(m) ,(m) 

Si - v^g - 1^5 


and 

s(>») = + rW = 1,2,3 (4.9) 

Where 

vi“\ ii'”^ - modal components of the fault-generated 

xs ■ xs . . 

components of the sending end phase voltages 
and line currents, determined by using the 
modal transformation matrix given in eqn, 
(4.7) • 
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«sot = R® ] 

and, 

= «iet = = Re[^3)] 

The amplitude comparison described in Section 4v2«l 
is applied to each signal pair, A tripping or blocking 
signal is issued by the first positive outcome of any one 
of the comparisons. Summarising, the tripping and blocking 
criteria are given by 

|c(®) ! ^|c(®) J A o f\ r^ .p ^ 14 . 

1^1 :• ^1^2 i»G*» ^ S^ ' - |S 2 I > 0,0 for faults 

, ahead 

and, 

<13^“^^ i.o., < 0,0 for faults 

behind, m = 1|,2,3 
(4.10) 

An identical relay arrangement is provided at the receiving 
end. 

4.2.3 Testing of the Relaying Principles by Digital Computer 
Simulation 

The viability of the new amplitude comparison techniques, 
described in earlier sections, is tested by considering the 
4-Bus system described in Chapter 3. In this system, all the 
lines are 128 Km in length. Frequency dependency of all the 
line and earth parameters is simulated in order to produce a 
realistic indication of the validity of the relaying principles. 
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An earth resistivity of ICX) ohm-m is assumed. Each 
equivalent source, representing the fault level at that bus, 
is assumed to have a ratio of zero sequence impedance to 
positive sequence impedance equal to 0,5 at power frequency, 
A phase sequence of a,b,c is assumed. 

The performance of the relay, assumed to be situated 
near bus 1 and affording protection to the line connected 
between buses 1 and 2, is studied for three-phase and one- 
phase- to-ground faults. For this purpose, the signal 
difference, viz,, | *"152^ t is computed for all the 

modes, and for the following locations of the fault, 

i) The remote end (marked F^) of the protected line? a 
distant internal fault 

ii) A point, 25,6 Km from the relaying point, on the 
protected line* a close-in internal fault 

iii) The remote end (marked F 2 ) of the line connected 
between buses 1 and 3} a distant external fault. 

In each case, fault resistances of 0,50,100,150 and 
200 ohms, and fault initiation angles of 0 and 90 degrees 
are considered. The signal difference for mode 1, in the 
case of 3-phase faults, has been found to be negligible. 
Therefore, ( - jS 2 |) modes 2 and 3 in the case of 

3-phase faults, and for all the modes in the case of one- 
phase— to— ground faults are plotted as a function of time. 
The plots are given in Figs, 4,3 to 4,32, 
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4*2.4 Discussion on the Results 

After a careful examination of the results, the 

following conclusions are drawn, 

i) The signal difference decreases with increase in the 
fault resistance. It depends also on the fault initia- 
tion angle. However, even for a. fault resistance of 
200 ohms and a fault initiation angle of 0°, it is 

of adequate value. The worst case is an external 
one-phase-to-ground fault, through R^ = 200 ohms 
and with 9^ = 0^, at the remote end of an adjoining 
line. In this case, the maximum signal difference is 
about -0,02 p,u, for mode 1* However, the signal 
difference of mode 2 is about -0,095 p,u, 

ii) For the 3-phase fault, only modes 2 and 3 are effective 
in producing tripping or blocking signals. This is as 
expected, since the mode 1 voltages and currents would 
be of negligible value by the very choice of the modal 
transformation matrix given in eqn, (4,7), 

iii) The results quite clearly confirm the viability of this 
simple relaying technique, 

4.3 FAULT LOCATING RELAY SCHEME 

This is the second UHS relay scheme,^ which is developed 

, ‘ ’I "I ' 

in this section, and which, fn addition to giving a tripping 

signal for internal faults, also gives the distance of the 
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fault. This scheme has the following special features. 

The tripping on internal faults takes place without the 
aid of communication channel. The communication channel 
is used only for blocking purposes in the event of external 
faults. Therefore, the operating time of the communication 
carrier equipments does not come into picture in the case 
of internal faults. The scheme makes use of the reflection 
properties of the travelling waves to determine the location 
of the fault. In this section^ the theory and principles of 
operation of the scheme, and also digital simulation results, 
that confirm the viability of the scheme, are presented, 

4,3,1 Theory and Principles of Operation 

With regard to this scheme also, the underlying 
principles are explained by considering a lossless 1-phase 
transmission line, which interconnects two large power 
systems, as shown in Fig, 4,1, Two identical relay schemes, 
located at S and R, and with reference directions as shown, 
are employed for the protection of the lino. The dis- 
cussion is centred on the relay scheme at S, It applies 
equally well to that at R. 

The input signal pair to the relay at S is the same 
as in the previous scheme, and is given by eqns, (4,3). 
However, in the present scheme, S 2 is used only to detect 
the reverse faults (i.e, faults behind the relay) whilst 

is used not only to detect but also to locate an internal 
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fault. These are explained below, 
(i) Reverse faults 


Upon the occurrence of a fault on any of the other 
lines connected to the sending end bus, travelling waves 
of voltage and current are generated. One pair of them 
travels towards the sending end bus where the waves undergo 
both reflection and refraction. Due to the refracted waves, 
|S 2 | increases from zero to Vf f ^ , as given 

by the second part of eqns. (4,4), where is the surge 
impedance of the faulted line. If ~ ^o’ usually 

the case, then | ^ff | ~ *^t '^ff | * There- 

fore, if I S 2 > 0, a reverse fault is indicated. However, 
there may be measurement and/or other types of errors. 
Accordingly, the following criterion is adopted for 
detecting a reverse fault. 


I S 2 I > a preset level (4,11) 

If j S 2 I exceeds the preset level, a transmitter is 
started after a short time delay Td^* The purpose of this 
time delay is explained later. In addition, issue of a trip 
signal at the local end, by the backward waves coming from 
the remote end after reflection, is blocked. The remote— end 
receiver, upon the receipt of the transmitter signal, blocks 
the issue of a trip signal there. 
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(ii) Forward faults 


Consider a fault at the point F (Fig. 4.i) on the 
protected line. The travelling waves, generated at the 
fault point, travel between this point and the ends of the 
line, for some time, due to multiple reflections and refra 
ctions. These events are depicted in Figs, 4.33 to 4.38, 
by means of Lattice diagrams, for different conditions. 


The first travelling wave pair reaching S will be 
backward and will cause the value of ^ increase 

suddenly from zero to ^}^f4 at this moment as given in 
Table 4,1, These waves will also cause | S2| to increase to 
Pg v^f I immediately thereafter. This may result in the 
transmitter being started. However, the startir^ of the 
transmitter is delayed by a short interval of time Td. , 

♦ X 

and meanwhile, the former event is used to block the 
starting of the transmitter. 


After a time Interval of 21:^^, there will be a second 
sudden increase in | S^| due to the backward reflected waves 
coming from the fault point. The time interval between the 
first two sudden increases in | | is measured. This time 

interval is egual to twice the travel time of the waves 
between *^0 fault and relay points. Hence, it gives the 
location of the fault. 


However, if the fault occurs at a voltage zero, there 


will not be any sudden increases in 


This fact is 
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FIG. 4-38 LATTICE DIAGRAM 

FOR A FAULT AT 

THE MIDDLE OF THE 
LINE 



2(1 + f|Pp)Vff 


max 


m ax 


Time 


F!G'4'39 VARIATION OF jS-jj WHEN . = 0 (not to scale) 
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illustrated in Fig, 4,39 where the time variation of 

shown for this case. This difficulty is surmounted 
by using sudden increase in d|s^|/dt in addition to those 
in| S^| for locating the faulty,: , i, , , , ,,,, 

For a close- in fault, | S^*j or djs^|/dt builds up 
rapidly and almost continuously. If the rapid and conti- 
nuous building up takes place over a certain preset time 
period, then tripping is initiated, and the fault is located 
as nearby* 


4*3.2 Relay settings 

In actual practice, the increases in either 




or 


d|Si|/dt will not be quite sudden, but will be only rapid. 
This is due to the attenuation caused by the losses in the 
line and earth. Also, the magnitude of increase in 
or d|s^|/dt depends upon the fault resistance which reduces 
'^ff* Accordingly, will have only a finite slope at 

the instant when the backward waves reach the relaying point, 
and not infinite. And, in between the first two rapid 
increases, |S^ I will have a slope equal to that of 
If =5 ^ffmax sin(«Qt+a), then the maximum slope of |S^| 
between these twi!Q,«t|istants could be 2 depending 

upon the fault initiation angle. The maximum value of 
will be fora bolted fault. Then, v^^ will be equal to 
( , the, profault voltage at the fault. And, the 

maximum slope of l^i/j would be equal to ^ ® o^pffjuax* 



165 


the rapid increase in | | is detected by measuring its 

slope, the following criterion needs adoption so that 
no wrong detection of fault occurs, when the slope of j S^| 
is maximum or close to it, during the time interval between 
the two rapid increases in|S^|, 

l.e. d|Sj^|/dt > (4.12a) 

At the instant of second rapid increase in|s^|, gets 
reduced to Hence, the following criterion is 

adopted for detecting the second rapid increase in | 




dt > 2 PgCii^ ^pffmax 


(4.12b) 


In a similar way, the following criteria are adopted 
to detect the first and second rapid increases in d|S^|/dt, 

d^Kj/dt^ > 2 0.2 


and 


■|s^|/dt" 


> 2 9 V 

s o pffmax 


(4.13) 


For very high fault resistances, the second rapid 
increase may not be appreciable* (This becomes evident when 
digital simulation results are observed). Under such cir— 

m''' 

cumstances, tripping may not ensue. This is avoided by the 
following procedure, A trip signal is issued after a fixed 
time delay of moment of detection of the 

first rapid increase either in | j or in d js^J/dt, provided 
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that no blocking signal is received from the other end. 

Of course, in such a case, fault location will not take 
place. 

The value of should be such that the following 
operations are not interfered with. In the event of an 
external fault beyond the remote end, the issue of a trip 
signal by the local relay should be blocked by the blocking 
signal received from the other end. For no interference with 
this operation, should be greater than the time taken 
by the communication equipment to produce the signal plus 
the response time of the blocking relay. In the event 
of an internal fault at the remote end of the protected 
line, the detection of second rapid increase in 
d|s^|/dt takes place after a time lapse, equal to twice the 
travel time of the entire line 2x, when reckoned from the 
moment of detection of first rapid increase. For no inter- 
ference with this operation, T ^2 ^ 

4.3.3 Typical cases 

Four typical cases of an internal fault are considered, 
in this section in order to illustrate the applicability of 
the basic principles explained in Section 4.3.1. 

(i) Bolted faults 

The Lattice diagram for this case is given in Fig. 4,33. 
It can be observed that |s^| increases suddenly from zero to 
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to 2 


V 


pff 


at T- , from 2 


lvfh° Vf| 


and so on» After sensing the first two sudden increases 
injs^l at and 3x^, a trip signal is issued. The time 
interval between these gives the location of the fault. 


The time of operation would be maximum for a fault at the 
remote end of the protected line, and would be equal to 3t. 
It obviously depends upon the length of the line. It would 
be equal to about 2 ms for a 200 Km line, and would be 
equal to 5 ms (i.e,, one-fourth of a cycle on 50 Hz 
frequency base) for a 500 Km line. 


(ii) Faults through resistance 

For this case, the Lattice diagram is given in Fig. 
4.34. An examination of the diagram reveals that, 
additional sudden changes in |S^| are produced by the 
transmitted (refracted) waves. The propagation of the 
transmitted waves is represented by broken lines on the 
diagram. However, these changes are normally decrements 
since P^, the voltage reflection coefficient at the fav^ 
point, is always negative, and and P^ are normally 
negative. Since only sudden increases in | | are sensed, 

these additional sudden changes do not affect the operation. 

(iii) Faults on an open-circuited line 

Let the receiving end be an open circuit. Then, the 
Lattice diagrams for (a) a bolted fault, (b) a fault through 
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closer to the sending end, and (c) a fault through 
closer to the receiving end would be as given in Figs. 4.35, 
4.36 and 4,37 respectively. It can be observed that normal 
operation takes place in the first two cases. In the third 
case, the sudden change at (t^+ 2 t: 2 ) would be an increment. 
Consequently, the relay operates faster. But, the fault 
point can bo located only at the second attempt as explained 
below. The fault is located first as at a distance of 
UT^ from the sending end, where u is the velocity of propa- 
gation of travelling waves. If the fault is not found there, 
then it is taken to be at a distance of 11 x 2 f^^om the receiv- 
ing end, 

(iv) Faults at the middle of the line 

For a fault at the middle of the line, there will not 
be a second or subsequent sudden increases in S^|, if 
p = P and Rv: = Z /2 as explained below, 

JL -*• w 

For this fault, the Lattice diagram will be as given 
in Fig. 4.38. When R^ = Zj2, 



From the Lattice diagram, it may be observed that, the 
second sudden increase at the sending end is given by, 

2[(1+ Pf) 



169 


With and = -0,5, the above quantity vanishes* 

That is, there will not be a second sudden increase. 

However, in the case of a 3-phase line, this will not 
pose any problem due to the following reason. The surge 
impedance of mode 1 and those of modes 2 and 3 have widely 
different values. And, if then or 

and vice versa. And hence, there will be at least 
one mode available for producing the desired operation, 

4,3,4 Extension of the Principles to the Protection of 
3-Phase Lines 

The extension of the relaying principles to the prote- 
ction of a 3-phase transmission line is based on the modal 
decomposition of the 3-phase line into three decoupled 
fictitious 1-phase lines, each having its own propagation 
constant and surge impedance, as discussed in Section 4,2,2, 
Consequently, the relaying principles described in Sections 
4.3.1 and 4,3,2 apply also to each of these fictitious 1-phase 
lines. Therefore, in this scheme also, three signal pairs, 
one for each mode, are employed. For the sending end relay, 
they are given by eqns, (4,9). The blocking and other 
criteria, given in eqns, (4,11) to (4,13), would get modified 
as given below in the same order. 

For reverse fault detection, 

> a preset level, m =1,2,3 (4,14) 
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For detecting 


.(m) 


(m 


the first and second rapid 
■ 1,2,3), we have 


increase in 


and 


.(m) 


/dt 


2(0 

o ffmax 




dTS 


1 / 


/(m) 


dt > 2Pg«Q ^Sax respectively; ra = 


1,2,3 

(4.15) 


For detecting the first and second rapid increases 
djs^’^^j/dt (m = 1,2,3), we have 


and 


(m) 

1 

(m) 


/dt® > 2 '^ 


ffmax 


/dt^ > 2Pg<.^V^®^^j^espectivelyi 


m 


in 


= 1,2,3 
(4.16) 


( m) 

Here, ® = 1,2,3 are the modal components of the 

fault generated components of the phase voltages at the fault 

point for a bolted fault. These components have different 

values for different types of faults as derived in the next 

section, and given in Table 4.3, From the table, it can be 

observed that, the value of is maximum in the case 

of one-phase-to-ground fault on phase ’a’ and is equal to 

( 2 ) ( 3 ) 

0,5 p.u. The values of ^fWax maximum for a 

3-phase fault, and each of them is equal to 0.57735G2 p.u. 

In order to provide safe settings, these maximum values are 
chosen for determining the settings defined in eqns, 

(4.14) and (4.15). V^f^ax ^ffraax will have lowest values 
for one-phasG-to-ground faults. Therefore, it is necessary 
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to verify that the settings chosen are. not detrimental in 
the case of one— phase— to-ground faults. The values of the 
settings for the detection of the first and second rapid 
increases are given in Table 4.2. 


TABLE 4.2 


Relay settings for the detection of first and second 
rapid increases in first and second time derivatives 

m = 1,2,3 


of S 


Nominal angular frequency of the system 

= 314.15926 rad/s 

fi = 0-^' 


Modes Detection of the Detection of the 

(m) first rapid increase in second rapid increase in 


first deriva- 
tive of 


second deri- first deri- 
vative of vative of 




second deri- 
vative of 
fgCm):, 


1 314.16 98696*04 251.33 78956.84 

2 and 3 362.76 113964.36 290.21 91171.49 


4.3.5 Modal Components of the Fault-Point Voltages for 
Different Types of Faults 

In this section, the modal components of the fault-point 
voltages are determined for different types of faults. 
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with the following assumptions*. 

1» The transmission line is assumed to be perfectly 
transposed 

2, The fault resistances are zero 

The conditions prevailing at the fault point are 
depicted in Fig. 4.40, where 
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(i) 3-phase fault 

For this fault, we have 


,a,b,c 

ff 


a , h j c __ . . 
pff ''pffmax 


Therefore, we have 


sin(tt^t+9^) 
s in ( t+(p^ -120® ) 
s i n (w ^ t+9^+ 120® ) 


1,2,3 _-l a,b,c 

^ff - ^ '^ff 


= V, 


pffmax 


(4.17) 


0.5773502 sin(«^t+<p^+30®) 
0.5773502 sin(wQt+9^-30®) 


(4.18) 


where S is the voltage modal transformation matrix, 
(ii) One-phase-to-ground-fault on phase a. 

For this fault, we have 


''If = ''Iff = pffmax 

and 

iff = iff =0 (4.19) 


Also, we know that (vide Appendix A) 

„1 , 2, 3 c*"^ * b,c 

'^ff - ^ '^ff 

and 

j1#2,3 

Iff - O Iff 


(4.20) 
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Solving eqns. (4,19) and (4,20), with S and Q having 
values as given in eqn. (4.7), we get 


•1,2,3 

^ff 


1 

1 

1 


i^ 


ff/3 


i.e., i^^^ = = 45^ (4.21) 




ff 


and [1 1 1] ^ ^ (4.22) 


Also, we have 


1%^*^ = + i5»^»^) 


ffs 


^ffr 


•1,2,3 

X • 0 • f X^^ — 

and = 


,/il»2,3 ^1,2,3\ 

^^ffS ^ ^-P-Pr* / 


^ffr 


yl,2,3 -1,2,3 _ ^1,2,3 .1,2,3 

o ^ffs “ o ^ffr 


(4.23) 

(4.24) 


where 


zJ-2>3 = diag[z(l) z'^^ Z<3)] 


For a perfectly transposed line, 
generally. 


( 2 ) _ 


and 


«5d. 2Z^^^ = 2Zp^ 

o o 0 


(4.25) 
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Solving eqns. (4,2l) to (4.25), we get 


V 


1,2,3 

ff 


V, 


pffmax 


0.50 sin((i»^t+9^) 
0.25 sin(t»)^t+q)^) 
0.25 sin(®Qt+g>^) 


(4.26) 


(iii) Phase-to-phase fault involving phases b and c 


For this fault, we have 


be be / b e \ 

''ff “ ''pff “ ^''pff ■ ''pff^ 


-V3 cos(.pt+9^) 


i^ — 0 
^ff ” ^ 


and 


= ° 


(4.27) 


Solving eqns. (4.27) in conjunction with eqns. (4.20) 
and (4.7), we get 


[1 - 1 ] 


V 


( 2 ) 

ff 


V 


(3) 

ff 


v’^^ V 

cos(cii 1+9 ) 

3 ^3 


id) 

Iff 


and [l l] 


i(2) 

Iff 

i(3) 

Iff 


= 0 


(4.28) 
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Finally, solving eqns. (4.28), (4.23), (4.24) and (4.25), 
we get 

(4.29) 

(iv) Double-phase- to-ground fault involving phases b and c 
For this fault, we have 

''ff = ''pff = sin("ot+To-120°) 

''ft = '^pff = ''ptt^ax 

and 

i|^ =0 (4.30) 

Solving eqns. (4.30), in conjunction with eqns. (4,20) 
and (4.7), we get 

[-2 1 1] ^ sin{<a„t+9^) 

and [1 1 1] = 0 (4.31) 

Finally, solving eqns. (4,31), in conjunction with eqns. 
(4,23) to ‘ (4.25) , we get 


,1,2,3 

'ff 


'pff max 


0 


0.2886751 cos(«Qt+9Q) 

-0.2886751 cos(“ t+m ) 

o 
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V 


1 , 2 , 

ff 


3 




pffmax 


. 0.400 sin(“^t+9^) 

0,3055 sin((i)^t+<p^-109.1®) 
0,3055 sin(w^t+9^4- 109.1*^) 


(4.32) 


The peak values of for the above faults, as 

obtained from eqns. (4.18), (4.26), (4.29) and (4.32), are 
tabulated in Table 4,3, 


Table 4.3 

Values of pffmax ~ lf2»3), for different types 
of bolted faults, expressed as a fraction of 


Type of fault 

''ffmax 

— 7 ^ 

ffmax 



ffmax 

Three-phase 

0,00 

0,5773502 

0.5773502 

0*ne~phase- to -ground 

0.50 

0,25 

0.25 

Phase-to-Phase 

0.00 

0,2886751 

0,2886751 

Do ubl e-phase- to - 
ground 

0,40 

0.3055 

0,3055 


4.3.6 Testing of the Relaying Principles by Digital 
Coraputer Simulation 

The viability of the relaying techniques, described in 
earlier sections, is tested by considering a 200 Km, 3— phase. 
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400 KV quad— conductor, single-circuit transmission line 
interconnecting two large power systems, similar to that 
shown in Fig, 4,1, The tower configuration and conductor 
spacings of the line are given in Fig, 3,7 and the basic 
data, in Appendix B, An earth resistivity of 20 ohm-m 
is assumed. The performance of the relay located near the 
sending end bus is studied for 3-phase and one— phase— to- 
ground faults. For simplicity, the following assumptions 
are made, 

i) The sending end and receiving end buses are fed 
exclusively by other lines and cables. Accordingly, 
the sending end and receiving end source networks are 
considered as purely resistive, 

ii) The resistances of the sending end and receiving end 
source networks are respectively equal to l/9th and 
l/4th of the real part of the appropriate modal surge 
impedance, 

iii) The transmission line is assumed to be perfectly 
transposed, 

iv) The parameters of the line are frequency invariant. 

The following values evaluated at the natural frequency 
of the line, viz,, 4,500 rad/s, are considered. 

Self and mutual series resistance/Km = 0,3992 and 


0.3476 ohms respectively 
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Self and mutual series indue tances/Km = 1,3215 and 

0,4198 mH respectively 

Self and mutual shunt capaci tances/Km =: 0,011299 and 

0,00187 |iF respectively, 

(v) The prefault peak voltage at the fault point is 1,0 p,u. 

By virtue of the assumption (iii) above, the relay 

(l) 

signal S^ would be zero for 3— phase faults and the signals 
(2) (3) 

S^ and Sj^ would be identical for one— phase— to— ground 
faults. Accordingly, the signals, S^^^ and for 3-phase 

faults, and and (is for one-phase— to— ground 

faults, are computed for the following locations of the 
faults, 

(i) At 50 Km from the sending end, 
and 

(ii) at the receiving end. 


In each of the above two cases, fault resistances of 
0,100,132,200 and 263 ohms and fault initiation armies of 0° 
and 90® are considered. The results are plotted in Figs, 

4,41 to 4,56, In Figs, 4,49, 4,50 , 4,53 and 4,54, the broken 
lines indicate the points where the slope of the curve has 
changed abruptly. 


Also, the values of d S^j/dt are computed for all the 
cases of the 3-phase fault, and for those cases of one-phase- 


to-ground faults having a fault initiation angle of 90°. The 
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First time derivatives of the relay signals in p.u. /second at the 
first rapid increase for a 3-phase fault 
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First time derivatives of the relay signals in p.u, /second 
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only in a very few cases, For the system studied, there is 
at least one unaffected relay signal, in these cases, to 
provide tripping. For example, in the case of a 3-phase fault 
at the receiving end and with m = 0°, and with = 200 
and 263 ohms, the second rapid changes in 


in j s‘3)^ 


are less 


than the settings. However, in these cases, the second 
rapid changes in are above the settings, and hence, 

the relay operates to give a trip signal. 


ii) The second rapid increases in the relay signals are 
low in the cases where 9^ =0®, and may fall below the relay 
settings for very high fault resistances. However, the 
occurrence of a fault at 9^ = 0® is relatively rare. In 
addition, the relay can still operate after a delay of T^2 
secQnds from the moment of detection of the first rapid 
increase, as explained in Section 4,3.2* Thus, this feature 
prevents the failure of the relay to operate under such 
relatively rare conditions. 


iii) The determination of the distance of the fault is 
illustrated with reference to Fig, 4,41, It can be observed 
that the time interval between the first two sudden increases 
in |s^| is about 325 ps. Therefore, the distance of the fault 
from the relay location is given by 

(3x10^) X (325x10“^) Ci* 48,75 km 
where 3x10^ is the velocity of the waves in km/s. 
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(iv) From the results given in Fig, 4,57, it can be 
observed that the relay would operate quite satisfactorily 
to give blocking signal for reverse faults. 

4.4 CONCLUSIONS 

In this chapter > two new relay schemes, based upon 
travelling wave phenomena, have been presented. The 
first scheme, which embodies the amplitude comparison of 
two relay signals, is simple and viable. The second scheme 
involves the use of the first or second time derivatives 
of a relay signal to determine the distance of the fault 
from the relaying point. In addition, the detection of an 
internal fault takes place without the aid of a communica- 
tion channel* This eliminates the time delays associated 
with the communication equipment. The determination of the 
distance of the fault, and the detection of internal faults 
without the necessity of communication between the ends of 
the protected line, are the two beneficial features of the 
second scheme. 

The viability of both the relay schemes, under various 
conditions, has been tested using digital simulation techniques, 
described in Chapter 3. The results obtained confirm their 
viability. 

The relaying schemes, based on travelling wave phenomena, 
are immune to power swings since power swings are comparatively 
slow transients, and in addition, these are seen as an external 
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disturbance. With these schemes, the detection of faults 
takes place before the current transformers (c.t.s) saturate 
and hence, the protection problems associated with the c.t. 
saturation are eliminated. Also, these relays operate 
before the protective gaps, across the series capacitors 
of the series-compensated transmission lines, flash over 
and thus, the protection problems associated with this 
occurrence are absent. 

The lightning strokes will have harmonic content in the 
J'jiHz range [72], Therefore, either the conventional trans- 
ducers with a cut-off frequency of a few KHz eliminate them, 
or a simple low-pass filter can be provided to eliminate 
them. The switching surges produced by switching operations, 
except those on the protected line, would be seen as produced 
by external faults. If the voltage transformers are placed 
on the line side of the circuit breakers, as has been the 
current practice on EHV systems, then the switching opera- 
tions, on the protected line also, are seen by the relays 
as external faults. In this case, the circuit breaker 
closure onto an existing fault on the protected line would 
be seen as an external fault. This difficulty can be over- 
come by a simple circuitry or logic as described in the 
reference [72], 



CHAPTER 5 


DIGITAL PROTECTION OF TRANSMISSION LINES USING TRAVELLING 

WAVE HiENOMENA 


5.1 INTRODUCTION 

With the emergence of inexpensive, fast and reliable mini- 
computers and microprocessors, the digital protection of trans- 
mission lines is going to be a reality in the near future. The 
digital protection of transmission lines provides an improved 
performance in terms of speed of operation as well as flexi- 
bility of obtaining, with ease, any desired relay characteri- 
stics. With suitable software logic, the protection schemes 
can be made self-checking against hardware failures. Also, the 
sequence of software events, which occurs in the processor in 
event of a fault, can be stored during the fault and output to 
a data link afterwards. 

A considerable amount of work in the field of digital 
protection of transmission lines has been reported since late 
1960’s. The development has been confined mainly to different 
types of algorithms which compute the fundamental frequency 
impedance. The algorithms proposed uptill now can be broadly 
classified into two groups j (l) distance (i.e. impedance) 
relay algorithms and (2) travelling wave relay algorithms. 
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The post-fault waveforms, during the first one or two 
cycles following the occurrence of a fault, comprise of a 
power frequency fundamental, an exponentially decaying d,c, 
component and high frequency transients. They also contain 
subtransient and transient power frequency components, if 
the line fault is near a generating source. The first group 
of algorithms involve in computing the power frequency 
impedance between the relaying and fault points from the 
fundamental components of the voltage and current obtained by 
suppressing other components through analog and/or digital 
filters and line modelling. The accuracy, with which the 
fundamental components of voltage and current can be deter- 
mined, depends upon the extent of filtering, and the extent 
and nature of filtering needed depends on the rigour of line 
modelling, both of which introduce longer time delays. Conse- 
quently, the speed and accuracy have been conflicting require- 
ments. On the other hand, the algorithms, based on travelling 
wave phenomena, make use of the complex post-fault waveforms 
with just one simple digital filtering to obtain the fault- 
generated components. Therefore, these algorithms would be 
accurate as well as fast. 

In this chapter, an overview of the distance relay algo- 
rithms is given first, in order to bring into focus their 
relative merits and inherent limitations. Next, two new 
travelling wave algorithms, whose viability has been tested on 



208 

a digital computer with the fault data generated by the 
digital simulation of sample power systems, are presented* 

5.2 AN OVERVIEW OF THE DIGITAL DISTANCE RELAY ALGORITHMS 

These algorithms involve in extracting the fundamental 
frequency components of voltages and currents from the 
complex post-fault waveforms and then, in determining the 
impedance between the fault and relaying points. Basically, 
there are four digital methods of determining the impedance 
from the fundamental components of voltage and current. In 
the first method, the magnitude of the impedance is computed 
as the ratio of the peak voltage to the peak current, and its 
argument, as the difference between the phase angles of the 
voltage and current. The peak values of, and the phase angle 
between, the voltage and current can be determined either from 
the samples of the fundamental components extracted from the 
ensemble of samples collected over one full power frequency 
period [88], or predicted from a much fewer number of samples 
of the fundamental components by using the samples and their 
time derivatives [89], or the first and second time derivatives 
of the samples [91] or the samples and the sampling interval 
[92]. 

In the second method, the fundamental components of vol- 
tage and current are determined in phasor form, from which the 
real and imaginary parts R and X of the pseudo impedance seen 
by the relay can easily be evaluated [87,93—100,107,111,112], 
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In the third method, the line is modelled by differential 
equations, the numerical solution of which yields the values 
of R and X [ 101-106, 113] • This method has the following 
advantages. It is necessary to filter out only those unwanted 
components v/hich are not covered by the line modelling. For 
example, if the line is modelled by a series R-L circuit, 
then there is no need to filter out the d.c. offset component. 
Also, this method permits the tackling, with ease, of the 
cases like the series compensated lines [113]. In the fourth 
method, the impedance is calculated directly in phasor form 
from the ratio of the frequency domain values of the fundamen- 
tal components of the voltage and current by the finite 
transform method [llO]. 

All those methods, which require samples over one full 
fundamental period for the computation of the pseudo impedance, 
do not offer any advantage in terms of speed when compared with 
the solid-state relays. The unwanted components of the post- 
fault waveforms are eliminated by employing various types of 
filters. Analog and digital filters, which have been proposed 
and/or used for this purpose, are described in the following 
sections. 

5,2.1 Analog Filters 

The mimic impedance used in the c.t. secondary plays the 
role of an analog filter since it filters out the d.c, offset 
in the current signal. Complete suppression of the d.c.offset 
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is, however, not possible owing to the fact that exact mat- 
ching of the X/R ratios of the primary and secondary circuits 
is difficult to obtain because the X/R ratio of the primary 
circuit upto the point of fault is variable. Some 
algorithms [89,91,97] assumed the use of this filter, 

RC low pass filters, with a suitable cutoff frequency, 
have been used to eliminate the high frequency components 
in conjunction with some of the algorithms [91,93], The 
fact, that a single analog filter cannot suppress all the 
unwanted components and that, these filters are slower than 
the digital ones, led to the development and preferential 
use of the digital filters, 

5,2,2 Digital Filters 

The various types of digital filters, which have been 
proposed so far, are described below. 

5. 2. 2.1 Notch Filters 

Two orthogonal notch filters with sine characteristics 
have been used to extract the fundamental components after 
the high frequency components are eliminated by an analog 
low pass filter [93], 

5 .2.2.2 Selected Harmonic Filter 

While the differential equations of the transmission line 
are being solved by numerical integration, the integration is 
carried out over a certain number of overlapping subintervals 
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with appropriate limits [lOSj. This leads to the elimination 
of certain harmonics and their multiples* 

5 *2,2, 3 Least Square Error Filters 

In one type, the deficiency in the modelling of the line 
by differential equations is treated as an error, and solu- 
tion for the parameters of the model is obtained subject to 
the minimisation of the mean of the squares of this error 
over the data window [102,103]. In another type, a polynomial 
fit (for example, a straight line fit over three points, a 
quadratic fit over five points or a cubic fit over seven 
points) is determined subject to the least square error 
criterion [89,97]. Differentiating this polynomial, the 
necessary time derivatives can be found. In yet another 
type, a waveform containing a decaying d.c. offset, the 
fundamental and a desired number of hamonic components is 
assumed, and the least square error criterion is applied to 
determine the unknown parameters of the fundamental component 
[107-109]. All these filters are rather slow and their 
accuracy depends on the data window as well as number of 
samples per cycle [ll9]. 

5. 2, 2, 4 Orthogonal Transform Filters 

Of these, the Fourier transform filter is the most 
widely used one, and utilizes sine and cosine functions 
as an orthogonal set [87,88,94-99]. This filter suppresses 
all the unwanted components, and therefore, offers the best 
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accuracy* But, it requires a data window; of one full 
fundamental period. However, filters employing half-a- 
cycle and less than half-a-cycle data windows with 
tolerable errors have been proposed [97,993* Filters using 
odd and even square waves [94], Walsh functions [lOOj, sample 
values and their derivatives [89j and the first and second 
derivatives of the samples [91] as orthogonal functions have 
also been proposed* 

5*2, 2, 5 Finite Transform Filter 

The fundamental frequency impedance is determined by 
carrying out the filtering process in -Uie frequency domain 
through the use of finite transfoim spectral analysis 
techniques [llOj* Acceptable accuracy limits could be 
achieved only with a speed of three-fourth of a cycle [llO]* 

5. 2. 2. 6 Kalman Filters [111,112] 

The nonfundamental components such as harmonics and d.c* 
components in the voltage and current waveforms are consi- 
dered as noise signals. The noise signal in the voltage 
waveform is considered as a white noise sequence with a 
decreasing variance while that in the current waveform as an 
exponential process plus a white noise sequence with a 
decreasing variance. Then, a two-state Kalman filter is used 
to extract the fundamental voltage phasor and a three- 
state Kalman filter to extract the fundamental current phasor. 
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The error in these filtering processes has been reported 
to be less than 1 percent after half-a-cycle [ill]. 

The filtering schemes, described so far, suppress the 
harmonics as well as d,c, offset from the post-fault current 
and voltage waveforms. However, in the case of ultra high 
speed relays, based on the travelling wave techniques, the 
transient components, such as high frequency harmonics and 
d,c, offset present in the post-fault current and voltage 
waveforms, are utilized for detecting the faults. As a 
consequence, these relays are faster and more reliable in 
operation. Hence, the filters, described above, will not be 
useful in the case of travelling wave relays. The cycle-to- 
cycle comparison scheme, proposed by Mann and Morrison [90 j, 
is adapted, in this work, as a digital filter for deriving 
the fault-generated components, which are used for forming 
suitable relay signals using the travelling wave principles* 

5.3 DIGITAL TRAVELLING WAVE RELAY ALGORITHf.1S 

This group of algorithms uses the fault-generated compo- 
nents of the relaying point phase voltages and line currents, 
obtained from the complex post-fault waveforms by means of a 
very simple digital filter and makes relayir^ decisions * 

through the application of travelling vmve techniques. This 
group offers two distinct advantages, which lead to a high 
speed of operation. First, there is no need for either 
analog or elaborate digital filters. Consequently, the time 
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delays associated with them are eliminated. Secondly, a 
very short data window can be employed. In fact, the 
sampling interval should be small enough to avoid the 
aliassing errors. That is, the sampling frequency should be 
atleast twice as great as the cutoff frequency of the 
transducers in order to avoid the aliassing errors. Thus, 
with a CVT of cutoff frequency of 2,5 KHz, the sampling 
frequency should be atleast 5 B3^z, and therefore, the 
sampling interval, 200 ps. Also, the number of samples 
needed to detect a fault would be small as explained later. 

With the high speed analog to digital conversion 
equipment and computers available at present, the above 
short intersampling period would be adequate to derive and 
store the six digital samples of the phase voltages and line 
currents, and to carry out the simple fault detection algo-» 
rithm described in the next section. 

Very few algorithms of this category have been proposed 
so far [71,114], most probably, due to the fact that the 
travelling wave relaying concept is relatively new. Two new 
algorithms of this category are presented in this chapter. One 
is based on the amplitude comparison relay and the other, on 
the fault locating relay, both described in Chapter 4, Each 
algorithm can be split into two parts, fault detection and 
relaying. The fault detection algorithm is common to both 
and, therefore, is described first. 
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5,4 FAULT DETECTION ALGORITHM 

This algorithm is similar to that proposed by Mann 
and Morrison [90], One counter is provided for each phase 
voltage. The difference between the currently sampled 
instantaneous voltage of each phase and the corresponding 
one in the previous cycle (which has been stored) is com- 
puted and stored. These differences, upto one cycle after 
the occurrence of a fault, would be the fault-generated 
components and are used in the relaying algorithms described 
later. If this difference is greater than a preset value, 
for any phase, then the counter of this phase is incremented 
by 1, OihervyisG, the concerned counter is decremented by 1, 
provided that it is not already zero. The latter ensures that 
no maloperation takes place due to spurious spikes. When 
the counter of any phase reaches a preset value, the fault 
is assumed to have occurred. Then, the sampling is sus- 
pended, the transmitter is switched off, and the relaying 
programme is executed. 

The setting, for the difference of instantaneous vol- 
tages one cycle apart, is chosen as 0.05 p.u. for both the 
relays. It is based on the assumption that load fluctua- 
tions do not lead to voltage variations of more. than 0.05 p.u. 
Actual value for this sotting can be determined by conducting 
digital Simula tio^i tests on the power system. The setting of 
each phase counter is taken as 5 for the amplitude comparison 



216 


relay, and as 10 for the fault locating relay. Higher 
setting for the latter is necessary to collect adequate 
number of samples. of fault-generated components of vol- 
tages needed for executing its relay algorithm. At any 
time, one full cycle of samples of voltages and currents 
are stored. As and when a new sample set arrives, the 
earliest set is deleted and the ensemble of sample sets is 
updated. 

It can be observed that this algorithm is quite simple, 
and therefore, can easily be carried out within an inter- 
sampling period. The flow chart of this algorithm is given 
in Fig, 5,1, 

5,5 AMPLITUDE COMPARISON RELAY SCHEME 

The theory and principles of operation of this scheme 
have been described in Chapter 4, The relaying -point phase 
voltages and line currents need sampling and then conversion 
to the digital form before they are supplied to the digital 
computer, in actual hardware implementation. From these, the 
fault-generated components are computed by using the cycle- 
to-cycle comparison method described by Mann and Morrison 
[90], and are used for relaying. The criteria, given in eqns, 
(4,10), are employed for detecting a fault ahead of and a 
fault behind the relaying point respectively. The computer 
application of this scheme, along with an algorithm, is des- 
cribed in the next section. 




KMTR(k) = KNTR(k)-1 


FALSE 
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5 •5*1 Digital Computer Application 

The modal components of the phase voltages and line 
currents are given by (vide Appendix A) 


and 






ja,b,c 


(5.1) 


As stated in Chapter 4, an ideally transposed transmission 
line is assumed for relaying purposes. With this assumption, 
S and Q become frequency independent and equal. And, 
choosing the Karrenbauer transformation matrix, given in 
eqn* (4fT), we can rewrite the above eqns. (5.1) as 


and 


where 


_ ^-1 ya,b,c 
f m f 

.1,2,3 _ .r“*l 4a,b,c 

f "" m f 

1 1 

1 -1 

1 0 




(5.2) 


(5.3) 


The computational effort is reduced by computing the modal 
components without division by 3, and these are called the 
modified modal components in the present work. Consequently, 
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the fault detection criteria get modified as given below, 

for faults ahead of the relaying 
point 

for faults behind the relaying 
pointj 

(5,4) 

where and are the modified relay inputs evaluated 

with the modified modal components of voltages and currents. 
With this modification incorporated, the relaying algorithm 
would be as follows. 


,,(m) 


,(m) 


and 


,,(m) 




m = 1,2,3 




Step 1 


Step 2 


Step 3 


Calculate, by the cycle-to-cycle comparison method, 
the fault-generated components of all the line 
currents for five sample sets backwards from the 
set at which the fault detection algorithm has 
yielded a logical 'yes* output. 

Calculate the modified modal components of the 
fault-generated components of the phase voltages 
and line currents. At the first instance, the 
earliest sample set is considered. 

Execute the Steps 4 to 6 for each mode, i,e,. 


m = 1,2,3 

Step 4 : Calculate the modified relay inputs 
Step 5 : Check if 

If so, a fault ahead of the relaying point has 
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has occurred. Then, start the transmitter* Also, 
if a carrier signal is received from the other end, 
trip the circuit breaker. If not, proceed to the 
next step. 

Step 6 : Check if 

If so, a fault behind the relaying point has 
occurred. Therefore, block the possible starting 
of the transmitter subsequently by the backward 
waves coming after reflection at the remote end of 
the line. If not, go to the next step. 

Step 7 : Take the next sample set and go to Step 2. 

The flow chart of the algorithm is given in Fig. 5.2. 

It can be seen that this algorithm is simple and straight- 
forward and therefore, can be implemented even with a micro- 
processor system. 

5.6 FAULT LOCATING RELAY SCHEME 

The theory and the principles of operation of this 
relaying scheme have been described in Chapter 4, The fault- 
generated components of the relaying-point phase voltages and 
lino currents are used for relaying purpose as in the previous 
scheme. VJith the modification . explained in the next section, 
the criterion given in eqn. (4.14) is used for detecting the 
reverse faults, and those given in eqns. (4.15) and (4.16) are 
used for detecting and locating internal faults. The computer 
application of this scheme, along with an algorithm, is des- 
cribed in the next section. 






!s 5, < So 


Tpjud Reverse fault has occurred 
^2 I y block the possible starting 
FALSE of the transmitter later ' 


200 ) 


No fault detected 
restrart sampling . 

FIG-5'2 flow chart OF THE AMPLITUDE 
COMPARISON RELAY ALGORITHM 
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5,6,1 Digital Computer Application 

The computation of division by 3, in evaluating the 
modal components of voltages and currents, is eliminated 
by using the modified components and the modified relay 
inputs, like in the previous scheme. Accordingly, the 
criterion for detecting the reverse faults is taken as 
given below, 

> 0.09 p,u,| m = 1,2,3 (5.5) 

And, the criteria for detecting the first and second rapid 


increases 

in 

for 

m = 1,2,3 are taken as. 


d /dt 

> 


and 

d s^(“^/dt 

> 

® " s“o fLx respectively, 


m = 1,2,3 (5.6) 


The criteria for detecting the first and second rapid 
increases in /dt for m = 1,2,3 get modified as, 

) 

and .. , 

j f""" , ' ' 

/dt^ > 6 Hfmax respectively? m = 1,2,3 

(5.7) 

I f 

, for 

m =1,2,3, at any sampling instant t^, are computed using 
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the following numerical differentiation formulae [120] 


d 


, (m) 


/dt] 


t=t . 


■Cl (ni) 

- jcf(m)i 

^l,i+l 

pi,i 


/At 


® “ 1»2,3 


and 

d2|s^(m) 




Cl (m) 

^l,i+l! 


= /(At)*", sayi m = 1,2,3 (5.9) 


-2 


±fX 


., (m) 

l,i-l 


(5,8) 


/ (At)^ 


where and are the modified relay 

inputs computed from (i+l)th, ith and (i-l)th sample sots of 
voltages and currents, and A t is tfte sampling interval. 

The first and, if necessary, the second derivatives are 
to be computed repetitively and hence, it is imperative that 
the computational effort be reduced. This is accomplished 
by modifying the criteria, used to detect the first and 
second rapid increases in , as follows* 

A > 6w At 

^1,1 o ffmax 


and 


^g^(m) At respectively} m = 1,2,3 

(5.10) 


Similarly, the criteria, used to detect the first and second 
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rapid increases in dfs 


i(ni) 


/ dt, are modified as given 


below« 


AC. (in) 

^ 1,2 




and 


A ct(in) 
^ 1,2 


> 6Pgft^ respectively! 


m = 1,2,3 

(5.11) 

As explained in Section 4,3,4, maximum values of 4?Lx 
m = 1,2,3 are chosen in the above settings. 


A rapid increase in Si 

? -*• 


(m) 


or 


:,(m) 


♦ 


/dt may span 


over two or three consecutive sampling intervals, in which 
case the fault-location computation may go wrong. This 
difficulty is overcome by a suitable logic, as described 
in the algorithm given below. With the above features 
incorporated j the relaying algorithm would be as follows. 


Stop 1 : Calculate, by the cycle-to-cycle comparison method, 
the fault-generated components of all the line 
currents for twelve sample sets backwards from 
the sample set at which the fault detection 
algorithm has yielded the logical ‘yes’ output. 

The two extra sample sets are used to partly cover the 
interval immediately preceding the arrival of waves 
at the relaying point in order to facilitate the 
computation of the second derivatives correctly. 
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Step 2 : Calculate the modified modal components of the 

fault-generated components of voltages and currents* 

Step 3 : Calculate modified relay inputs for all the modes 
and for all the sample sets. 

Step 4 ; Check if > 0,09 p.u. for m = 1,2,3. 

If so for any m, stop executing the relaying 
programme further at this end and start the 
transmitter, after a time delay of Td^ (defined 
in Section 4.3,1 of Chapter 4), which sends infor- 
mation to the other end to stop executing the 
relaying programme there also and to block tripping. 
If not, proceed further. At the start, this step 
is executed for the first sample set. 

Step 5 i Calculate ASl^?^ and check if . 

Xt-L ■ XfX O 

^ffmax “ 1>2,3, If not, take the next 

set of modified inputs and go to the Step 4. If 
satisfied, do not go to Step 4 any further. Block 
the possible starting of the transmitter and set 
the time counter to the value of the time corres- 
ponding to this sample set. The time counter 
sotting and the rapid-increase-detection criterion 
remain the same if the above condition is satisfied 
over consecutive sampling intervals. If this 
happens for four* consecutive sampling sets, an 
indication, that the fault is close-in, is given as 



226 


output, otherwise, change the setting of the rapict- 

increase~detection criterion to 6p u At. 

•^s o ffmax 

When the changed criterion is satisfied for any m, 
indicating a second rapid increase, set the time 
counter to the current time minus the previous setting. 
Stop further execution of the algorithm, compute the 
distance of the fault and output this value as well as 
a logical 'yes’ for tripping. When this step yields 
no result, go to the next step. 

Step 6 : Repeat the algorithmic procedure of Step 5 with the 

second derivatives of .S^^’'”^; for m = 1,2,3, using the 
rapid-increasG-detection criteria given in eqns. (5.11). 
If this step also fails to yield any result, go to the 
next step. 

Step 7 : Under these conditions, the second rapid increase 

detection must have failed. Then, a trip signal is 
issued after a time delay of Td 2 » which is defined in 
Section 4,3,2 of Chapter 4. 

The flow chart of the algorithm is given in Fig. 5,3. 
5.7 DIGITAL SIMULATION AND RESULTS 
5,7,1 Testing 

For testing the proposed algorithms, fault data are 
generated for a desired type of fault for two cycles. 
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one prefault and one post-fault, by employing the digital 
simulation techniques described in Chapter- 3* The 
amplitude comparison relay algorithms are tested by using 
the fault data computed for an internal fault at the 
remote end (marked F^) of the protected line and for an 
external fault (marked F 2 ) at the remote end of the 
adjoining lino in the 4-bus system shown in Fig. 3.11. 

Each source is represented by an equivalent lumped- 
parameter model with phase inductances calculated from the 
fault level specified and the neutral impedance, from an 
assumed ratio of ZjZ^ of 0.5 at power frequency, as 
explained in Appendix D. A three-phase fault through a 
resistarice of 100 ohms per phase and with a fault initia- 
tion angle of 0® is assumed. 

The fault locating relay algorithms are tested by using 
the fault data evaluated for a single line to ground fault 
on phase ’a’ at the remote end of the protected line of 
200 Km length in the power system shown in Fig, 4,1, In 
this case, the sending - and receiving-end sources are 
considered as composed solely of transmission lines and 
cables, and as having net modal surge impedances respectively 
equal to l/9th and l/4th of the appropriate modal surge 
impedances of the line, A fault resistance of lOO ohms per 
phase and a fault initiation angle of 90^ are considered. 
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In both the sample power systems, each 3-phase trans- 
mission line is a typical 400 KV quad-conductor untransposed 
single-circuit line. An earth resistivity of 100 ohm-m 
is assumed and the frequency dependency of all the line 
and earth parameters is taken into consideration. The 
remaining data of each line are given in Fig. 3.7 and 
Appendix B. The sampling interval is taken as 200 ps and 
the peak value of the prefault voltage at the point of 
fault is assumed as 1.0 p.u. 

Using the realistic fault data, generated as described 
above, the fault detection and relaying programmes are run 
on the DEC System 1090 computer at I.I.T,, Kanpur. 

5.7.2 Results 

In both the cases, the CPU time has been found to be less 
than 1 ms for the execution of the entire programme, 
consisting of fault detection and relaying. However, eight 
sample sets in the first scheme and fourteen in the second 
were required for detecting the internal fault. Hence, the 
total time of operation would be about 8x0. 2+1 = 2.6 ms in 
the first scheme, and it would be about 14x0.2+1 = 3.8 ms in 
the second scheme. The second scheme takes more time 
primarily on account of the fault locating feature. Also, 
in the second scheme, the distance of the fault has been 
computed as 210.0 Km, the actual value being 200 Km, Thus, 
this scheme can locate the faults only approximately. 
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5*7,3 Discussion 

Summarising, the digital travelling wave algorithms, 
developed in this chapter, have been tested in the following 
manner* The fault data are generated for a sample power 
system by employing realistic models for the transmission 
lines, and the accurate frequency-domain method. Using 
these data as input, the programmes of the relay algorithms 
are executed on the digital computer, DEC System 1090 at the 
I.I.T*, Kanpur, for testing the viability of the above 
algorithms. Since the object has been only to test the 
viability of the algorithms, only one faulted condition is 
considered in each case. 

However, in actual hardware implementation, that is, 
in the on-line digital computer application, the three phase 
fault voltages and currents, derived from the system through 
transducers, have to be sampled simultaneously at the 
desired sampling frequency, converted to the digital form, 
and transferred to a digital processor, by means of an 
analog subsystem. The currents are monitored as voltage drops 
across resistors on the secondaries of the interposing current 
transformers. A typical analog subsystem, for this purpose, 
may consist of six buffer amplifiers, six sample— and— hold 
circuits, a multiplexer and an analog-to-digital converter 
(ADC). The purpose of the buffer amplifiers is to trim the 
input signals. The digital output from the ADC is transferred 



231 


to the digital processor, that is, the computer, through a 
direct memory access (DMA) channel. All the components of 
the analog subsystem would have to be of high-speed type* 

A simplified block schematic diagram of the hardware, needed 
ill actual implementation, is given in Fig. 5.4. The hardware 
implementation, which would yield more realistic indication 
of the viability of the computer application of the travelling 
wave relay algorithms, could not be carried out for want of 
facilities. 

5.8 CONCLUSIONS 

The digital simulation tests, carried out on the two 
new digital travelling wave relay algorithms described in 
this chapter, have established their viability. The heart of 
these schemes is the fault detection algorithm, which is quite 
simple and straightforward. With this fault detection 
algorithm, not only the relaying algorithms described in this 
chapter but also the analog travellirig wave relay schemes 
proposed by others [73,77,85] can be implemented digitally. 
There being no need for any elaborate filtering, the travelling 
wave relay algorithms are superior in terms of both the speed 
and accuracy. Therefore, these will ultimately supplant 
the distance algorithms in the digital protection of trans- 
mission lines. 



Buffer amplifier modules 















CHAPTER 6 


CONCLUSIONS 


6.1 GENERAL 

The continual growth, in the size and complexity 
of electric power systems, has given rise to the need 
for extremely fast and highly reliable protective schemes. 

The utilization of the travelling wave phenomena has enabled 
the realization of u.h.s. relaying schemes for the protection 
of EHV/UHV transmission lines. And, the application of 
digital computers for protective relaying purposes is 
replete with several advantages, j the chief amongst l^em 
being the self-checking facilities and the maximum flexi- 
bility. Therefore, the use of travelling wave phenomena, 
along with a digital computer, for the protection of 
EHV/UHV transmission lines, has potential merits. Accord- 
ingly, the primary objective of this thesis has been the 
development of new relaying schemes based on the travelling 
wave phenomena, and also, the exploration of the feasibility 
of the adaptation of these schemes for the digital computer 
application. In the following sections, a brief account of 
the work carried out in this thesis, and the scope for 
further work, are presented, 

6.2 REVIEW OF THE WORK CARRIED GUT IN THIS THESIS 

Use of solid-state circuitry renders the generation 
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of several special relay threshold characteristics, like 
conic and quadrilateral etc*, easier and simpler, and 
also the attainment of high speed of operation possible. 

And, the use of ICs, in place of transistors and several 
discrete components, results in more compact, more 
reliable and economical relay units. The quadrilateral 
pickup characteristic is eminently suitable for the 
protection of long and heavily-loaded EHV/UHV transmission 
lines. Therefore, new distance relays, employing ICs 
and capable of generating a three-step quadrilateral pickup 
characteristic, have been developed, fabricated and tested. 
The first relay scheme demonstrates that the use of ICs 
greatly simplifies the circuitry, although the input 
signals used are the same as those given in reference [ 36 ], 
The second relaying scheme uses less number of input 
signals for generating the open quadrilateral characteristic, 
and therefore, is more reliable. Also, the measurement of 
the coincidence of the input signals is completed in half 
a cycle, and there are no time co-ordination problems 
associated with this scheme. Results of the tests 
conducted on these relays closely agree with the theoretical 
ones. 

For designing, and for assessing the performance of, 
both analog u.h.s. relay schemes and high-speed digital 
computer relay algorithms, especially those based on 



236 


neodod for not more than one cycle* Under such circum- 
stances, the frequency-domain method of computation does 
not need undue computational effort. Also, the application 
of the fast Fourier transform method to the computation of 
time-domain values does not provide much computational 
relief. The results of numerical examples, worked out in 
this thesis, establish these facts* 

Ultra high speed clearing of faults improves the 
transient stability* The fault clearing time depends on 
the speed of the protective relay as well as on that of 
the associated circuit breaker* With the emergence of 
u,h*s, circuit breakers, the need for u*h,s. protective 
relays has become imperative* The development of these 
relays has been facilitated by the use of travelling wave 
phenomena. However, only a fev/ relays, based on travelling 
wave phenomena, have been developed so far* Consequently, 
new travel ling -wave relay schemes, suitable for the 
protection of EHV/UHV transmission lines, have been 
developed in the present work. 

In one relaying scheme, the amplitude comparison of 
each of the three pairs of modal relay inputs is made use 
of to distinguish between the reverse and forward faults. 
Tripping is initiated if the fault is detected as forward 
at both the ends of the line. The underlying relaying 
principles have been verified by the results obtained 
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through the digital simulation of a sample power system* 

In the second scheme, one relay input (for each iMde) 
is used for detecting the reverse faults, and tripping is 
blocked at the remote end with the aid of a carrier comiau-* 
nication channel and also at the local end, under reverse 
fault conditions* By making use of the first, and if 
necessary the second, time derivatives of the other relay 
input (for each mode), the instants of first and second 
incidences of the backward travelling waves, at the relaying 
point, can be determined* This would enable the determina- 
tion of the distance of the fault and hence, the detection 
of an internal fault* Thus, both the detection of an 
internal fault and the determination of the distance of the 
fault can be accomplished with this relaying scheme* Since 
carrier communication between the ends of the line is not 
necessary for the detection of internal faults, the operating 
time of the carrier comraunication equipment does not form a 
part of the operating time of the relaying scheme, unlike 
in the other travelling-wave relaying schemes. The 
applicability of the fault-locating principle has been 
illustrated by considering four typical cases of an 
internal fault. In addition, the relaying principles 
involved in this scheme have been validated by the results 
of digital simulation studies conducted on an equivalent 
sample power system. 
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The effect of fault resistances and fault initiation 
angle on the relay response, in both the schemes, has been 
studied* The response of both the relaying schemes has 
been found to be good, even for fault resistances as high 
as 200 ohms. However, the response may be poor for higher 
fault resistances when the fault initiation angle is zero. 
However, the occurrence of a fault at a voltage zero is 
comparatively rare* With the application of relaying 
principles to each of the three propagation inodes, all 
the types of faults are taken care of in both the schemes. 
Therefore, these relaying schemes can be regarded as a novel 
type of polyphase relays. 

The development of the digital protection of transmission 
lines has been confined mainly to the generation and improve- 
ment of different types of distance relay algorithms, which 
involve in the computation of the fundamental power frequency 
impedance from the complex post-fault current and voltage 
waveforms. Different types of analog and digital filters 
have been developed for extracting the fundamental components 
of voltages and currents. But, the accuracy, with which the 
fundamental components can be determined, depends upon the 
extent of filtering. On the other hand, the nature and 
extent of filtering needed depends upon the rigour of line 
modelling. Either more elaborate filtering or more rigorous 
line modelling would introduce longer time delays. As a 
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consequence* the speed and accuracy have been conflicting 
requirements with the digital distance algorithms. But, 
the digital algorithms, based on travelling wave phenomena, 
necessitate the filtering out of only the prefault compo- 
nents. Also, a fewer number of samples at shorter inter- 
vals of time are required for fault detection. Thus, 
these algorithms provide faster relaying. 

In this thesis, the cycle-to-cycle comparison method, 
proposed by Mann and Morrison [90], has been adapted as a 
simple digital filter for filtering out the prefault compo- 
nents. Also, algorithms, for the digital computer application 
of the amplitude-comparison and fault-locating relay schemes, 
have been developed and presented. The viability of these 
algorithms has been tested on the digital computer, DEC 
System 10^ at I.I.T., Kanpur, by using realistic fault 
data obtained from the digital simulation of sample power 
systems. The results of these tests have demonstrated the 
viability of the proposed algorithms. However, there is a 
need for on-linedigital computer testing of these algorithms, 
Vi/hich could not be carried out for want of facilities, to 
establish, more firmly, the viability of these algorithms. 

6.3 SCOPE FOR FURTHER WORK 

The fault data, for testing the u.h.s, protective relays 
developed in this thesis, have been obtained by the digital 
simulation of sample power systems, wherein, although the 
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transmission lines have been modelled with adequate details 
and accuracy, the sources considered are of a simple nature. 
Therefore, the frequency-domain representation of complex 
source networks incorporating generation is clearly an area 
in which further work needs to be done. In addition, the 
frequency-domain models of current and voltage transducers 
as well as power transformers need incorporation in the 
system model to include their transient behaviour. 

These will facilitate in the generation of more realistic 
fault data. 

The validity of the travelling-wave relaying principles, 
developed in this thesis, needs further assessment v\4\en applied 
to the protection of double-circuit lines, series and shunt 
compensated long lines and multi-terminal lines. 

Actual hardware implementation of, and extensive field 
tests on, the travelling-wave relay schemes, proposed in this 
thesis, are necessary to warrant their satisfactory perfor- 
mance under actual conditions. 


I 



APPENDIX A 


FREQUENCY DOMAIN A,B,C,D PARAJ4ETERS OF A 
TRANSMISSION LINE 

A three-phase transmission line with distributed 
parameters, is described in the frequency domain by the 
following differential equations [116] • 

dx^ 

and, = IF.b.ojt^a.b.c 

dx^ 

where 

V^'bfC column vector of the transforms of the voltages 
to earth of the phases a,b,c at a point X on the 
line, 

■ja,b,c ^ column vector of the transforms of the line 

currents, pertaining to the phases, a,b,c, at X 
p,b,c ^ r^.b.C, :^,b,C 

Z^'bjC^ t^,b,c ^ transforms of the p*u, length series impedance 

and shunt admittance matrices, respectively, of 
the line 

and, X — distance of X from S*E. 
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and are evaluated by taking into account 

the presence of earth wires, and the frequency dependency of 
both the resistive and inductive parameters of the line and 
earth* The theory of natural modes enables the above 
equations to be effectively decoupled by means of 
of voltage and current modal transformation matrices S and Q 
respectively. Using these transformation matrices, equations 
(A,l) can be rewritten as follows. 


^251,2.3 ^ 


T7 

dx 


and 


dx 


(A.2) 


where 


Id o 

S - the matrix formed from the eigenvectors of P ' * 

Q = , (i,e* the matrix formed from the ; eigen- 
vectors of ) 

y 1»2,3 _ g-1 ^,b,c^ column vector of modal components of 

phase voltage transforms 


= q”^ column vector of modal components of 

line current transforms 
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S-l [pa.b.C] s = Q-1 [^.b.cjt ^ ^ ^ 

and and X^ - eigenvalues of 

The equations (A,2) can be rewritten in the decoupled 
form as given below. 


= » vf">) = v2 v^”) 

I m 


dx 


m 


and. 


^ -1(111) ^ 2 -(m) „ = 1 2 ■ 

dx 


'm 


(A.3) 


These equations show that the wave propagation in a 
3-phase transmission line can be decomposed into three 
separate and independent components, called modal components, 
each possessing its own propagation constant and surge 
impedance i 


Solving the first equation of (A,3), we get 

Y^'^) _ cosh Yjh^ + ^2 sinh m = 1,2,3 (A#4) 

where and K 2 are arbitrary constants to be determined by 
boundary conditions. 

Also, it is known [II 6 ] that 

^a , b , c 7a,b,c 

‘h ^ 


^a,b,c 

dx 
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Upon transformation to the modal domain, the above equation 
becomes 


^1.2,3 ^ ^1,2,3 

v/h ere 

\ = diag. Z^^^] 


(A.5) 


The equation (A*5) can be rewritten in the decoupled 
form as given below. 


dx “ ^ 


m = 1,2,3 


(A.6) 


Solving equations (A, 4) and (A.6), we get 

= -Ifer sinh Yiu^ + ^ cosh ^ = 1*2,3 

^ (A.7) 

Now, we know that, at x = 0, i.e., at the S.E,, 
and at x = 1, i,e., at the R.E., 

and ■= 


Using the above boundary conditions and solving eqns. (A,4) 
and (A, 7) , we get 


Km) 


= (cosh Y^l) - sp (sinh Y„l). 


and 

T(m) 


■R - - ffm) ■" ■= 1'2,3 


t(">), .r, „ 


(A. 8) 
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Expressed in a matrix form, equations (A, 8) become 


■^1,2,3 

R 

1 


[coshyl] 

-^[Yj^HsinhYl] 


yl f 3 

: 

■^R 


“[yJ^^CsinhYlJ 

[coshyl] 


t1,2,3 

2 


(A. 9) 

where 


[coshyl] 

= diag,[coshY^l 

coshY2l 

coshy^l] 

[sinhyl] 

= diag,[sinhy2^1 

sinhY2l 

sinhy^l] 

Cy] 

= diag. [y^ y^ 

^3^ 



From.eqn, (A. 9), the following equation can easily be 
derived. 


“^1,2,3 

• 

[coshyl] 

^}yiCY]”^[ sinhyl] 


^1,2,3 

Y1,2,3 

J-s _ 


[yJ ^^[ sinhyl] 

[coshyl] 




*^, 2,3 


» say 


(A.IO) 
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Performing inverse modal transformation on the above 
equation, vve get 





APPENDIX B 


BASIC DATA OF THE SINGLE-CIRCUIT 400 KV 3-PHASE 
QUAD-CONDUCTOR TRANSMISSION LINES 


Number of conductors per phase 

Number of earth wires 

Conductor resistivity 

Earth-wire resistivity 

Conductor strand diameter 

Earth-wire strand diameter 

Geometric mean radius of the 
bundle conductor 

Number of effective strands 
in phase conductors 

Number of effective strands in 
the earth wire 

Outer diameter of earth wire 


= 4 
= 1 

= 3.20x10"® Ohm-m 
- 2.69x10”® Ohm-m 
= 0.32 cm 
= 0.32 cm 

= 15.495 cm 

= 54 

= 54 

= 2.86 cm 



APPENDIX C 


SOURCE DATA OF THE 3-BUS SAMPLE POWER SYSTEM 


Nominal angular frequency = 314*159 rad/s 
Phase sequence = a,b,c 
Source impedances : 

r,_ = r,.u = r_. = O 

c 

L. 

K 

L. 

L. 


aa 

= ^bb 

= He = 

aa 

^bb 

II 

O 

O 

it 

ab 

” He 

^ca 

’ba 

JQ 

O 

it 

~ He 



APPENDIX D 


DERIVATION OF EQUIVALENT SOURCE MODELS FROM SPECIFIED 

FAULT LEVELS 


A simplified equivalent-source model, based upon 
the given short-circuit level at a bus, would be as shown in 
Fig. D.l. The phase impedances are calculated from the fault 
level specified, and the neutral impedance, from a specified 
ratio of zero and positive sequence source impedances as 
outlined briefly below [ll5]. 


Let the fault level at a bus be x MVA, 
the line voltage of the bus be V KVj 
the ratio of the zero and positive sequence 
source impedances be 

the nominal frequency of the system be rad/s. 


Then, the phase impedances would be given by 
Z = = Z = Z_ 


V^ . 

-a -b - “c - “ph = - 


(D.l) 


It can easily be shown that, 
Z„^+3Z 
'ph 




3Zn 
n ^ 

^ph 


{D.2) 


where Zp^^ and are the phase and neutral impedances res 
pectively. 




FIG. D-1 SIMPLIFIED EQUIVALENT SOURCE 
MODEL 
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Solving equations (D,l) and (D,2), we get 

Zn = - ^ ohms (D.3) 

Neglecting resistances, we get 

% 

Sh “ 7 ^ 

and 

Ln = - ^ ^ 

Equations (D.4) and (D..5) have been used in connection 
with the 4-bus sample power system, employed in the thesis. 
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